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SUMMARY 


A  theoretical  research  program  directed  toward  the  study  of  the  energe¬ 
tics  and  LWIR  radiative  properties  of  selected  uranium/oxygen  band  systems  has 
been  undertaken.  Included  in  this  research  program  was  the  investigation  of 
the  strongest  electronic  and  vibrational  bands  in  the  LWIR  region  for  the 
species  UO,  U0+,  UC^ ,  U02  +  ,  and  UC>2++.  The  program  for  accomplishing  this 
research  effort  was  formulated  into  three  separate  tasks:  a)  adaption  of  our 
electronic  structure  codes  to  the  DNA  CYBER  76  System,  b)  calculation  of  per¬ 
tinent  electronic  wavefunct ions  and  energies,  as  a  function  of  internuclear 
separation  and  within  a  relativistic  framework,  for  selected  species  of  the 
uranium/oxygen  system  which  may  be  important  in  the  LWIR  region,  and  c)  calcu¬ 
lation  of  electronic  transition  moments  and  transition  probabilities  between 
specific  vibrational  levels  of  the  electronic  states  corresponding  to  the 
strongest  radiating  band  systems  belonging  to  the  uranium/oxygen  system  and 
prediction  of  IR  and  possible  optical  oscillator  strengths. 

Our  calculations  indicate  that  the  species  U0+  will  be  efficiently  solar 
pumped  and  will  exhibit  strong  radiation  in  the  region  0.6  <  X  <  11.3  p. 
Further,  we  predict  efficient  conversion  of  solar  photons  to  IR  photons  for 
this  species.  For  U02  +  ,  very  strong  solar  pumping  is  predicted  for  X  <  0.5  p. 
The  pathways  for  conversion  of  solar  photons  to  the  LWIR  region  are  still 

uncertain  but  strong  LWIR  radiation  (fQ1  ~  10-4)  is  predicted  for  the 
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vibrational  transitions  of  the  ground  state  of  UO2  .  Further  studies  of 

UO  and  examination  of  the  relative  importance  of  the  doubly  ionized  species 
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(U  ,  UO  ,  UO2  )  in  this  system  are  indicated. 


PREFACE 


This  report  was  prepared  by  the  United  Technologies  Research  Center,  East 
Hartford,  Connecticut,  under  Contract  DNA  001-82-C-0015.  The  research  was 
performed  under  Program  Element  62715H,  Project  125AAXH,  Task  Area  x,  Work 
Unit  00043  and  was  funded  by  the  Defense  Nuclear  Agency  (DNA). 

Inclusive  dates  of  research  were  1982  January  6  through  1983  March  31. 

Mr.  Peter  W.  Lunn  (RAAE)  was  the  Contract  Technical  Manager  (CTM)  for  this 
contract . 

Very  useful  discussions  with  Drs.  Russell  Armstrong  (Air  Force  Geophysics 
Lab),  Morris  Krauss  (National  Bureau  of  Standards),  Robert  Field  (MIT)  and 
G.  A.  Peterson  (UTRC)  are  also  acknowledged. 

All  aspects  of  the  research  work  reported  herein  were  aided  by  the 
skilled  help  of  Judith  B.  Addison  (UTRC),  who  carried  out  much  of  the  computer 
program  development  and  assisted  in  the  analysis  of  the  calculated  data  and  in 
the  preparation  of  this  final  report. 


TABLE  OF  CONTENTS 


Section  Page 

SUMMARY  .  1 

PREFACE  .  2 

LIST  OF  TABLES .  4 

1  INTRODUCTION  .  5 

2  METHOD  OF  APPROACH  -  NONRELATIVISTIC  METHODS  .  7 

A.  Quantum  Mechanical  Calculations  .  7 

1.  Levels  of  Approximation  .  7 

2.  Spin  and  Symmetry .  7 

3.  Method  of  Ab  Initio  Calculation  .  8 

4.  Molecular  Integrals  .  11 

5.  Configuration  Selection  .  11 

6.  Density  Functional  Approach  -  Method  ...  12 

7.  Computational  Aspects  of  the  XQ  Method  ....  16 

B.  Transition  Probabilities  .  20 

3  DISCUSSION  OF  RELATIVISTIC  METHODS  .  30 

A.  Breit-Pauli  Hamiltonian  .  32 

B.  Approximate  Treatments  .  33 

C.  Effective  Core  Models .  37 

4  DISCUSSION  OF  RESULTS  .  41 

UO .  41 

U0+ .  43 

U02 .  45 

-U02+ .  46 

U02++ .  47 

5  RECOMMENDATIONS  .  48 

REFERENCES .  73 


LIST  OF  TABLES 


Molecular  states  of  UO  . 

Calculated  spectroscopic  data  for  UO/UO+  .  .  . 
Calculated  oscillator  strengths  (f^,^,)  for  the 
vibrational  -  rotational  transition  of  UO  (X5I 
Total  integrated  absorption  coefficients  for 

UO  ground  state  . 

Molecular  states  of  U0+  . 

Term  levels  for  U+1(5f37s2),  cm-1 . 

Calculated  oscillator  strengths  (f^,^,,)  for  the 
vibrational  -  rotational  transition  of  U0+  (A4I 
Total  integrated  absorption  coefficients  for  UO' 
ground  state  . 


SECTION  1 


INTRODUCTION 

The  release  of  certain  chemical  species  into  the  upper  atmosphere  results 
in  luminous  clouds  that  display  the  resonance  electronic-vibration-rotation 
spectra  of  the  chemically  reacting  species.  Such  spectra  are  seen  in  rocket 
releases  of  chemicals  for  upper  atmospheric  studies,  upon  re-entry  into  the 
atmosphere  of  artificial  satellites  and  missiles,  and  as  a  result  of  energy 
deposition  in  the  atmosphere  caused  by  nuclear  weapons  effects.  Of  particular 
interest  in  this  connection  is  the  observed  spectra  of  certain  metallic 
oxides.  From  band  intensity  distributions  of  the  spectra,  and  knowledge  of 
the  f-values  for  electronic  and  vibrational  transitions,  the  local  conditions 
of  the  atmosphere  can  be  determined  (Reference  1).  Such  data  are  fundamental 
for  the  analysis  of  detection  and  discrimination  problems. 

Present  theoretical  efforts,  which  are  directed  toward  a  more  complete 
and  realistic  analysis  of  the  transport  equations  governing  atmospheric 
relaxation  and  the  propagation  of  artificial  disturbances,  require  detailed 
information  of  thermal  opacities  and  LWIR  absorption  in  region  of  temperature 
and  pressure  where  both  atomic  and  molecular  effects  are  important  (References 
2  and  3).  Although  various  experimental  techniques  have  been  employed  for 
both  atomic  and  molecular  systems,  theoretical  studies  have  been  largely 
confined  to  an  analysis  of  the  properties  (bound-bound,  bound-free  and  free- 
free)  of  atomic  systems  (References  4  and  5).  This  has  been  due  in  large  part 
to  the  unavailability  of  reliable  wavefunctions  for  diatomic  molecular 
systems,  and  particularly  for  excited  states  or  states  of  open-shell 
structure.  Only  recently  (References  6-8)  have  reliable  procedures  been 
prescribed  for  such  systems  which  have  resulted  in  the  development  of 
practical  computational  programs. 

The  application  of  these  computational  methods  to  studies  of  the  elec¬ 
tronic  structure  and  radiation  characteristics  of  metal  oxides  has  been 
reported  for  several  of  the  lighter  systems  (References  9-11).  A  preliminary 
study  of  the  uranium/oxygen  system  has  been  reported  by  Michels  (Reference 
12)  which  identified  a  large  number  of  low-lying  molecular  states  for  both  the 


UO  and  U0+  systems.  Of  particular  interest  was  the  discovery  of  two  struc¬ 
tures  for  U0+  that  resulted  from  two  different  spin-couplings  of  the  uranium 
valence  electrons.  These  results  suggested  strong  LWIR  radiation  in  U0+  aris¬ 
ing  from  different  electronic  state  transitions. 

An  inherent  uncertainty  in  these  preliminary  calculations  was  present, 
owing  to  the  neglect  of  relativistic  effects  that  were  much  too  difficult  to 
include  in  molecular  calculations  at  that  period  of  time.  The  7s  valence 
electron  of  uranium,  and  its  corresponding  o-bonding  molecular  orbital,  are 
highly  relativistic  in  nature  which  results  in  a  contracted  charge  density 
relative  to  that  which  would  occur  in  lighter  molecular  systems.  The  effect 
of  this  contraction  on  the  relative  positions  of  the  low-lying  electronic 
states  of  the  uranium/oxygen  system  can  now  be  calculated  with  some  degree  of 
confidence  using  newly  developed  relativistic  computer  codes. 

Because  of  inherent  difficulties  in  the  experimental  determination  of  the 
spectroscopy,  transition  probabilities  and  LWIR  radiation  for  metal  oxide 
systems  and  in  light  of  the  aforementioned  recent  progress  in  the  calculation 
of  relativistic  electronic  wavefunctions ,  especially  for  diatomic  systems,  a 
technical  program  for  calculating  these  properties  was  undertaken  for  the 
Defense  Nuclear  Agency  under  Contract  DNA001-82-C-0015.  The  emphasis  in  this 
work  was  on  the  ions  of  uranium  and  uranium  oxides  (U  ,  UO  ,  UC^)  since  these 
have  been  determined  to  be  important  radiators  in  the  LWIR  region.  These 
studies  indicated  that  the  doubly  ionized  species,  U+  ,  U0++  and  UO**,  should 
also  be  considered  because  of  their  role  in  charge  neutralization  processes 
and  their  potential  as  early-time  radiators. 

The  general  composition  of  this  report  is  as  follows.  In  Section  II,  we 
present  a  description  of  the  mathematical  methods  which  were  employed  in  this 
research.  Included  in  Section  II  are  sub-sections  which  deal  with  the 
construction  of  electronic  wavefunctions,  the  calculations  of  expectation 
properties,  the  evaluation  of  molecular  transition  probabilities,  and  the 
calculation  of  electronic  wavefunctions  using  both  the  ab  initio  and  density 
functional  methods.  This  is  followed  by  Section  III  which  describes  the 
inclusion  of  relativistic  effects  into  the  density  functional  (Xa)  method. 

The  calculated  results  and  pertinent  discussion  are  presented  in  Section  IV. 
Recommendations  are  presented  in  Section  V. 


SECTION  2 


METHOD  OF  APPROACH  -  NONRELAT1VISTIC  METHODS 
A.  QUANTUM  MECHANICAL  CALCULATIONS 

Central  to  these  theoretical  studies  are  the  actual  quantum-mechanical 
calculations  which  must  be  carried  out  for  the  atomic  and  molecular  species. 
For  added  clarity,  various  aspects  of  these  calculations  are  discussed  in 
individual  subsections. 

1.  Levels  of  Approximation 

Much  evidence  on  diatomic  and  polyatomic  systems  indicates  the  inadequacy 
of  a  minimum  Slater-type-orbital  (STO)  basis  for  constructing  quantitatively 
correct  molecular  wavefunctions  (References  13  and  14).  This  means  inner- 
shell  and  valence-shell  STO's  of  quantum  numbers  appropriate  to  the  atoms  (Is, 
2s,  2p,  for  C,  N,  0;  etc.).  The  main  deficiency  of  the  minimum  basis  set  is 
in  its  inability  to  properly  describe  polarization  and  the  change  of  orbital 
shape  for  systems  which  exhibit  large  charge  transfer  effects.  Values  of  the 
screening  parameters  4  for  each  orbital  can  either  be  set  from  atomic  studies 
or  optimized  in  the  molecule;  the  latter  approach  is  indicated  for  studies  of 
higher  precision.  When  high  chemical  accuracy  is  required,  as  for  the 
detailed  studies  of  the  ground  or  a  particular  excited  state  of  a  system,  a 
more  extended  basis  must  be  used.  Double-zeta  plus  polarization  functions  or 
optimized  MO's  are  required  for  reliable  calculated  results  of  chemical 
accuracy. 

The  chosen  basis  sets  give  good  results  only  when  used  in  a  maximally 
flexible  manner.  This  implies  the  construction  of  Cl  wavefunctions  with  all 
kinds  of  possible  orbital  occupancies,  so  that  the  correlation  of  electrons 
into  overall  states  can  adjust  to  an  optimum  form  at  each  geometrical 
conformation  and  for  each  state.  Except  when  well-defined  pairings  exist,  as 
for  closed  shell  and  exchange  dominated  systems,  a  single-configuration  study 
(even  of  Hartree-Fock  quality)  will  be  inadequate. 

2.  Spin  and  Symmetry 

Proper  electronic  states  for  systems  composed  of  light  atoms  should 
possess  definite  eigenvalues  of  the  spin  operator  S2  as  well  as  an  appropriate 
geometrical  symmetry.  The  geometrical  symmetry  can  be  controlled  by  the 
assignment  of  orbitals  to  each  configuration,  but  the  spin  state  must  be 


obtained  by  a  constructive  or  projective  technique.  Formulas  have  been 
developed  (Reference  15)  for  projected  construction  of  spin  states  from 
orthogonal  orbitals,  and  programs  implementing  these  formulas  have  been  in 
routine  use  at  UTRC  for  several  years. 

One  of  the  least  widely  appreciated  aspects  of  the  spin-projection  prob¬ 
lem  is  that  the  same  set  of  occupied  spatial  orbitals  can  sometimes  be  coupled 
to  give  more  than  one  overall  state  of  given  S  quantum  number.  It  is  neces¬ 
sary  to  include  in  calculations  all  such  spin  couplings,  as  the  optimum  coup¬ 
ling  will  continuously  change  with  changes  in  the  molecular  conformation. 

This  is  especially  important  in  describing  degenerate  or  near-degenerate 
excited  electronic  states. 

3.  Method  of  Ab  Initio  Calculation 

A  spin-free,  nonrelativistic ,  electrostatic  Hamiltonian  is  employed  in 
the  Born-Oppenheimer  approximation.  In  systems  containing  atoms  as  heavy  as 
Kr,  this  approximation  is  quite  good  for  low-lying  molecular  states.  For  a 
diatomic  molecule  containing  n  electrons,  the  approximation  leads  to  an  elec¬ 
trostatic  Hamiltonian  depending  parametrically  on  the  internuclear  separation, 

<^(R)=-f  i -^-i ^  ♦  ♦  z  t-  (1) 

2  i*  I  '  i*l  'A  i-.i  '  i  B  R  i>j>,  'I  KU 

where  and  Zg  are  the  charges  of  nuclei  A  and  B,  and  r^A  is  the  separation 
of  electron  i  and  nucleus  A.  is  in  atomic  units  (energy  in  hartrees,  length 
in  bohrs) . 

Electronic  wavefunctions  j»(R)  are  made  to  be  optimum  approximations  to 
solutions,  for  a  given  R,  of  the  Schrodinger  equation 

^(R)t(R)=  E  (R)i(R)  (2) 

by  invoking  the  variational  principle 


8W  (R)  =  8 
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The  integrations  in  Equation  (3)  are  over  all  electronic  coordinates  and 
the  stationary  values  of  W(R)  are  approximations  to  the  energies  of  states 
described  by  the  corresponding  States  of  a  particular  symmetry  are 

studied  by  restricting  the  electronic  wavefunction  to  be  a  projection  of  the 
appropriate  angular  momentum  and  spin  operators.  Excited  electronic  states 
corresponding  to  a  particular  symmetry  are  handled  by  construction  of  config¬ 
uration-interaction  wavefunct ions  of  appropriate  size  and  form. 

The  specific  form  for  iJi(R)  may  be  written 

*  (R)  =  Z  (R)  (4) 

H- 

where  each  ’t'^(R)  is  referred  to  as  a  configuration,  and  has  the  general  struc¬ 
ture 

=  n(V  (5) 


where  each  is  a  spatial  orbital,  /)  is  the  antisymmetrizing  operator, 


is  the  spin-projection  operator  for  spin  quantum  number  S,  and  9  is  a  product 

M 


of  a  and  3  one-electron  spin  functions  of  magnetic  quantum  number  M.  No 
requirement  is  imposed  as  to  the  double  occupancy  of  the  spatial  orbital,  so 
Equations  (4)  and  (5)  can  describe  a  completely  general  wavefunction. 


In  Hartree-Fock  calculations  ip(R)  is  restricted  to  a  single  ^  which  is 


assumed  to  consist  as  nearly  as  possible  of  doubly-occupied  orbitals.  The 


orbitals  are  then  selected  to  be  the  linear  combinations  of  basis  orbitals 


best  satisfying  Equation  (3).  Writing 


$u.  i  -  £  a  v  i  X  i 


(6) 


the  aV£  are  determined  by  solving  the  matrix  Hartree-Fock  equations 


x.o.i  (each  X) 


(7) 


where  e-  is  the  orbital  energy  of  <J>  ;. 


The  Fock  operator  has  been  thoroughly  discussed  in  the  literature 
(Reference  16)  and  depends  upon  one-  and  two-electron  molecular  integrals  and 
upon  the  ay^.  This  makes  Equation  (7)  nonlinear  and  it  is  therefore  solved 
iteratively.  UTRC  has  developed  programs  for  solving  Equation  (7)  for  both 
closed  and  open-shell  systems,  using  basis  sets  consisting  of  Slater-type 
atomic  orbitals.  Examples  of  their  use  are  in  the  literature  (Reference 
7). 

In  configuration  interaction  calculations ,  the  summation  of  Equation  (4) 
has  more  than  one  term,  and  the  c^  are  determined  by  imposing  Equation  (3)  to 
obtain  the  secular  equation 


where 


Z  ws^„)c„  =  o  (each/u.) 

HM1/  =/ (R )J  (R)¥„(R)d/r 
S (R)V*  (R)tfT 


Equation  (8)  is  solved  by  matrix  diagonalizat ion  using  either  a  modified 
Givens  method  (Reference  17)  or  a  method  due  to  Shavitt  (Reference  18). 

The  matrix  elements  and  may  be  reduced  by  appropriate  operator 

algebra  to  the  forms 

=  Z  «p  ^sP  'fy'  «^(R)P  n  ¥„j  (Xi.R^  (10) 

v  .  \  /n  n  V  (11) 

€P \®m  p  ®M/\n (i I  ,R)  P  n  *„j  (Xi  ,rz> 

where  P  is  a  permutation  and  its  parity.  The  sum  is  over  all  permutations 
<0MkVly  *s  3  "Sanikel  coefficient"  and  the  remaining  factors  are  spatial 
integrals  which  can  be  factored  into  one-  and  two-electron  integrals.  If  the 
are  orthonormal,  Equations  (10)  and  (ll)  become  more  tractable  and  the 
H^v  and  S,,„  may  be  evaluated  by  explicit  methods  given  in  the  literature 


(Reference  15).  Computer  programs  have  been  developed  for  carrying  out  this 
procedure,  and  they  have  been  used  for  problems  containing  up  to  106  total 
electrons,  10  unpaired  electrons,  and  several  thousand  configurations. 

The  Cl  studies  described  above  can  be  carried  out  for  any  orthonormal  set 
of  <(>^  for  which  the  molecular  integrals  can  be  calculated.  Programs  devel¬ 
oped  by  UTRC  make  specific  provision  for  the  choice  of  the  $  .  as  Slater-type 
atomic  orbitals,  as  symmetry  molecular  orbitals,  as  Hartree-Fock  orbitals,  or 
as  more  arbitrary  combinations  of  atomic  orbitals. 

4.  Molecular  Integrals 

The  one-  and  two-electron  integrals  needed  for  the  above  described  method 
of  calculation  are  evaluated  for  ST0*s  by  methods  developed  by  the  present 
investigators  (Reference  19).  All  needed  computer  programs  have  been  devel¬ 
oped  and  fully  tested  at  UTRC. 

5.  Configuration  Selection 

Using  a  minimum  basis  plus  polarization  set  of  one-electron  functions,  a 
typical  system  can  have  of  the  order  of  104  configurations  in  full  Cl  (that 
resulting  from  all  possible  orbital  occupancies).  It  is  therefore  essential 
to  identify  and  use  the  configurations  describing  the  significant  part  of  the 
wavefunction.  There  are  several  ways  to  accomplish  this  objective.  First, 
one  may  screen  atomic-orbital  occupancies  to  eliminate  configurations  with 
excessive  numbers  of  anti-bonding  orbitals.  A  third  possibility  is  to  carry 
out  an  initial  screening  of  configurations,  rejecting  those  whose  diagonal 
energies  and  interaction  matrix  elements  do  not  satisfy  significance  criteria. 
Programs  to  sort  configurations  on  all  the  above  criteria  are  available  at 
UTRC. 

Other,  potentially  more  elegant  methods  of  configuration  choice  involve 
formal  approaches  based  on  natural-orbital  (Reference  20)  or  multiconfigura¬ 
tion  SCF  (Reference  21)  concepts.  To  implement  the  natural-orbital  approach, 
an  initial  limited-CI  wavefunction  is  transformed  to  natural-orbital  form,  and 
the  resulting  natural  orbitals  are  used  to  form  a  new  Cl.  The  hoped-for 
result  is  a  concentration  of  the  bulk  of  the  Cl  wavefunction  into  a  smaller 
number  of  significant  terms.  The  multiconfiguration  SCF  approach  is  more 
cumbersome,  but  in  principle  more  effective.  It  yields  the  optimum  orbital 
choice  for  a  preselected  set  of  configurations.  This  approach  works  well  when 
a  small  number  of  dominant  configurations  can  be  readily  identified. 
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It  should  be  emphasized  that  the  probLem  of  configuration  choice  is  not 
trivial,  and  represents  an  area  of  detailed  study  in  this  research.  The  exis¬ 
tence  of  this  problem  causes  integral  evaluation  to  be  far  from  a  unique 
limiting  factor  in  the  work. 

6.  Density  Function  Approach  -  Xa  Model 

The  Xa  model  (Reference  22)  for  the  electronic  structure  of  atoms,  mole¬ 
cules,  clusters  and  solids  is  a  local  potential  model  obtained  by  making  a 
simple  approximation  to  the  exchange  -  correlation  energy.  If  we  assume  a 
nonrelativistic  Hamiltonian  with  only  electrostatic  interactions,  it  can  be 
shown  that  the  total  energy  E  of  a  system  can  be  written  exactly  (Reference 
23)  (in  atomic  units)  as 


E  =  Z  Dj  <  Uj 


•  /* 


U.  > 


+  4-1 


*u*v 


fJL^U  TfJ.U 
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This  expression  is  exact  provided  the  uj  are  natural  orbitals  and  n^  are 
their  occupation  numbers  (i.e.,  eigenfunctions  and  eigenvalues  of  the  first 
order  density  matrix).  The  first  term  in  Equation  (12)  represents  the  kinetic 
and  electron-nuclear  energies.  The  second  term  is  the  nuclear  repulsion 
energy.  The  sums  (p,v)  are  over  all  the  nuclear  charges  in  the  system.  The 
third  term  is  the  electron-electron  repulsion  term,  which  represents  the  clas¬ 
sical  electrostatic  energy  of  the  charge  density  p  interacting  with  itself, 
where 


^>(l)  =Z  rvu^tO  UjO) 


(13) 


The  last  term  Exc  represents  the  exchange  correlation  energy  and  can  be 
expressed  formally  as 

Exc  =~rfp{i)  d'T /  (l;2)  drT  ,  (14) 

12 

where  Pxc  (1,  2)  represents  the  exchange-correlation  hole  around  an  electron 
at  position  1.  In  the  exact  expression,  pxc  is  dependent  on  the 


12 


is  the 


second-order  density  matrix.  In  the  Hartree-Fock  approximation  E 
exchange  energy,  p  represents  the  Fermi  hole  due  to  the  exclusion  principle 
and  depends  only  on  the  first-order  density  matrix.  In  the  Xa  method,  we  make 
a  simpler  assumption  about  p  .  If  we  assume  that  the  exchange-correlation 
hole  is  centered  on  the  electron  and  is  spherically  symmetric,  it  can  be  shown 
that  the  exchange-correlation  potential 


Uxc  =  / 


drT 

r  12  2 


(15) 


is  inversely  proportional  to  the  range  of  the  hole,  rs,  where  rs  is 
defined  by 


47 r 
3 


s  P 


(I)  =  I 


(16) 


Therefore,  in  the  Xa  model,  the  potential  is  proportional  to  p^^(r).  We 

define  a  scaling  parameter  a  such  that 


UH  (')=  -  nr"  (3/>(l)/87r)'/3 


(17) 


The  expression  in  Equation  (37)  is  defined  so  that  a  =  2/3  for  the  case  of  a 
free  electron  gas  in  the  Hartree-Fock  model  (Reference  24)  and  a  =  1  for  the 
potential  originally  suggested  by  Slater  (Reference  25).  A  convenient  way  to 
choose  this  parameter  for  molecular  and  solid  state  applications  is  to  opti¬ 
mize  the  solutions  to  the  Xa  equations  in  the  atomic  limit.  Schwarz  (Refer¬ 
ence  26)  has  done  this  for  atoms  from  z  =  1  to  z  =  41  and  found  values  between 
2/3  and  1. 

In  the  "spin  polarized"  version  of  the  Xa  theory,  it  is  assumed  (as  in 
the  spin-unrestricted  Hartree-Fock  model)  that  electrons  interact  only  with  a 
potential  determined  by  the  charge  density  of  the  same  spin.  In  this  case  the 
contribution  to  the  total  energy  is  summed  over  the  two  spins,  s  =  ±  1/2. 


where  the  potential  is  spin-dependent 


and  ps  is  the  charge  density  corresponding  to  electrons  of  spin  s.  The 
spin  polarized  Xa  model  is  useful  for  describing  atoms  and  molecules  with 
open-shell  configurations  and  crystals  which  are  ferromagnetic  or  anti¬ 
ferromagnetic  . 

Once  one  has  made  the  Xa  approximation  to  the  total  energy  functional  E 
in  Equation  (12),  then  the  rest  of  the  theory  follows  from  the  application  of 
the  variational  principle.  The  orbitals  u£  are  determined  by  demanding  that 
E  be  stationary  with  respect  to  variations  in  u£.  This  leads  to  the  set  of 
one-electron  Xa  equations 


u.  *  «j  Uj 


(20) 


where  is  the  one-electron  eigenvalue  associated  with  u^  .  Since  p(r)  is 
defined  in  terms  of  the  orbitals  u^,  Equation  (20)  must  be  solved  itera¬ 
tively,  until  self-consistency  is  achieved.  Empirically,  if  one  takes  as  an 
initial  guess  that  p  is  approximately  a  sum  of  superimposed  atomic  charge 
densities,  then  the  convergence  of  this  procedure  is  fairly  rapid.  The  factor 
of  2/3  multiplying  the  potential  is  a  result  of  the  linear  dependence  of  E^c 
on  p.  This  also  has  a  consequence  that  the  Xa  eigenvalues  e£  do  not  satisfy 
Koopman's  theorem,  i.e.,  they  cannot  be  interpreted  as  ionization  energies. 
However,  it  can  be  shown  that  the  are  partial  derivatives  of  the  total 
expression  of  Equation  (12)  with  respect  to  the  occupation  number, 


If  E  were  a  linear  function  of  n£,  then  Koopmans'  theorem  would  hold. 

However,  because  of  the  dominant  Coulomb  term,  E  is  better  approximated  by  a 
quadratic  function  in  n£.  This  leads  to  the  "transition  state"  approxima¬ 
tion  which  allows  one  to  equate  the  difference  in  total  energy  between  the 

state  (n. ,  n.)  and  (n.-  1,  n.+l)  to  the  difference  in  the  one-electron 
i  J  i  J 

energies  e.  -  e.  calculated  in  the  state  (n.  -  1/2,  n.  +  1/2).  The  error  in 
3  1  i  J 

this  approximation  is  proportional  to  third-order  derivatives  of  E  with 

respect  to  n£  and  nj ,  which  are  usually  small  (Reference  27).  The  main 

advantage  of  using  the  transition  state  rather  than  directly  comparing  the 

total  energy  values  is  computational  convenience,  especially  if  the  total 

energies  are  large  numbers  and  the  difference  is  small. 

The  relationship  of  Equation  (21)  also  implies  the  existence  of  a  "Fermi 

level"  for  the  ground  state.  This  can  be  seen  by  varying  E  with  respect  to 

n-  under  the  condition  that  the  sum  Z  n^  is  a  constant,  i.e., 


(e-xI 


(22) 


implies  3E/ 3n£  ■  X,  where  X  is  a  Lagrangian  multiplier.  This  implies  that 
the  total  energy  is  stationary  when  all  the  one-electron  energies  are  equal. 
However,  the  occupation  numbers  are  also  subject  to  the  restriction  0  <  n£  <  1 
This  leads  to  the  following  conditions  on  the  ground  state  occupation 
numbers ; 

«.  <X )  n.  =  I 

i  i 

«i  >  X  )  rig  =0  (23) 

€ j  =  X  ■  )■  0<  r»j  <  i 

In  other  words,  the  ground  state  eigenvalues  obey  Fermi  statistics  with  X 
representing  the  Fermi  energy.  It  should  be  noted  that,  in  contrast  to  the 
Hartree-Fock  theory,  where  all  the  n£  are  either  0  or  1 ,  the  Xot  model  pre¬ 
dicts,  in  some  cases,  fractional  occupation  numbers  at  the  Fermi  level.  In 
particular,  this  will  occur  in  a  system  (such  as  transition  metal  or  actinide 
atom)  which  has  more  than  one  open  shell. 


The  Xa  model  differs  in  other  significant  ways  from  the  Hartree-Fock 
method.  In  fact,  the  simplification  introduced  in  approximating  the  total 
energy  expression  introduces  several  distinct  advantages  over  Hartree- 
Fock: 

1.  The  primary  advantage  is  purely  computational.  The  one-electron 
potential  in  Equation  (20)  is  orbital-independent  and  local,  i.e.,  it  is  the 
same  for  all  electrons  (except  in  the  spin-polarized  Xot  theory)  and  is  a 
multiplicative  operator.  On  the  other  hand,  the  Hartree-Fock  potential  is 
nonlocal,  or  equivalently,  there  is  a  different  local  potential  for  each 
orbital.  This  involves  a  great  deal  more  computational  effort,  especially  for 
systems  described  by  a  large  number  of  orbitals.  It  has  been  shown  (Reference 
28)  that  the  Xa  orbitals  for  the  first  and  second  row  atoms  are  at  least  as 
accurate  as  a  double-zeta  basis  set,  and  are  probably  better  for  larger  atoms 
which  involve  electrons  with  Z  >  2. 

2.  The  orbital-independent  Xa  potential  leads  to  a  better  one-electron 
description  of  electronic  excitations  of  a  system.  Both  the  unoccupied 

(n^  =  0)  and  occupied  (n^  =  1)  eigenfunctions  are  under  the  influence  of  the 
same  potential  resulting  from  the  other  N-l  electrons.  The  Hartree-Fock  vir¬ 
tual  orbitals  see  a  potential  characteristic  of  the  N  occupied  orbitals,  and 
therefore  are  not  as  suitable  for  describing  the  excited  states.  Actually, 
although  the  ground  state  virtual  eigenvalues  are  usually  a  good  description 
of  the  one-electron  excitations,  the  virtual  spectrum  of  the  transition  state 
potential  where  one-half  an  electron  has  been  removed  from  the  system  gives  a 
much  better  first-order  picture  of  these  levels  (Reference  29). 

3.  As  has  been  shown  by  Slater  (Reference  30),  the  Xa  model  rigorously 
satisfies  both  the  virtual  and  Hel Iman-Feynman  theorems,  independent  of  the 
value  of  the  parameter  a.  This  is  convenient  for  calculating  the  force  on  a 
nucleus  directly  in  terms  of  a  three-dimensional  integral,  rather  than  the 
six-dimensional  integrals  in  the  expression  for  the  total  energy  of  Equation 
(12). 

7.  Computational  Aspects  of  the  Xa  Method 

In  application  of  the  Xa  model  to  finite  molecular  systems,  there  are  two 
practical  aspects  of  the  calculations  which  must  be  considered.  The  first 


concerns  the  choice  of  the  integration  framework  for  describing  the  molecular 
wavefunctions  and  the  second  deals  with  the  choice  of  the  exchange  parameter, 
o,  in  different  regions  of  space. 

In  computations  with  heteronuclear  molecules,  there  are  several  free 
parameters  that  must  be  chosen:  the  ratio  of  sphere  radii  for  the  atomic 
spheres  of  integration  at  a  given  internuclear  separation,  the  degree  of 
sphere  overlap,  and  the  value  of  the  exchange  parameter  in  the  atomic  spheres 
and  the  intersphere  region. 

It  has  been  found  that  changing  the  ratio  of  the  sphere  radii  for  the  two 
atoms  in  a  heteronuclear  diatomic  molecule  introduces  changes  in  the  total 
energy  that  can  be  large  on  a  chemical  scale  (~  1  eV) .  A  choice  for  sphere 
radii  based  on  covalent  bonding  radii  does  not  necessarily  provide  a  good 
estimate  for  these  calculations.  The  value  of  the  exchange  parameter,  a,  and 
the  sphere  radii  and/or  sphere  overlap  is  normally  fixed  in  Xa  calculations 
for  crystals  where  the  geometry  is  fixed.  However,  to  develop  a  potential 
curve,  the  molecule  description  needs  to  change  substantially  as  the  inter¬ 
nuclear  separation  varies  and  the  changing  sphere  radii  include  varying  frac¬ 
tions  of  the  total  molecular  charge  (Reference  31).  Studies  made  at  UTRC  have 
shown  that  at  any  given  separation  the  total  energy  calculated  from  the  Xa 
model  is  a  minimum  at  the  radii  ratio  where  the  spherically  averaged  poten¬ 
tials  from  the  two  atomic  centers  is  equal  at  the  sphere  radius. 


V1(rg.)  =  V9(r, 
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(24) 


This  relationship  between  the  potential  match  at  the  sphere  boundary  and  the 
minimum  in  the  total  energy  appears  to  hold  exactly  for  "neutral"  atoms  and 
holds  well  for  ionic  molecular  constituents.  In  the  case  of  two  ionic 
species,  the  long  range  tail  of  the  potential  must  go  like  +2/R  from  one  ion 
and  -2/R  (in  Rydbergs)  for  the  other  ion  and  so  at  large  internuclear  separa¬ 
tions,  the  tails  of  the  potential  cannot  match  well.  However,  at  reasonable 
separations,  the  1/R  character  of  the  potential  does  not  invalidate  the  poten¬ 
tial  match  criterion  for  radii  selection.  This  match  for  the  atomic 
potentials  is  applied  to  the  self-consistent  potentials. 


In  molecules  with  significant  charge  sharing  in  the  bonds,  the  radii  of 
the  atomic  spheres  is  frequently  increased  in  Xa  calculations  so  that  an  over¬ 
lap  region  appears  in  the  vicinity  of  the  bond  (Reference  32).  Studies  made 
at  UTRC  show  that  the  contribution  to  the  total  molecular  energy  from  the 
exchange  integral  shows  a  minimum  at  the  optimum  sphere  radius  or  sphere  over¬ 
lap.  This  provides  a  sensitive  criterion  for  selecting  these  parameters. 

The  values  of  the  exchange  parameters  in  the  spherical  integration  region 
around  each  atomic  center  are  frequently  set  at  the  atomic  values  both  for 
neutral  and  for  ionic  molecular  constituents.  However,  for  light  atoms,  the 
value  of  a  which  best  reproduces  Hart ree-Fock  results  varies  substantially 
with  ionicity.  In  argon,  the  following  table  compares,  for  the  neutral  atom 
and  the  positive  ion,  the  HF  energy  and  the  Xa  energy  calculated  for  several 
values  of  a. 


a 

Xa  Energy 

HF  Energy 

Ar° 

.72177 

526.8176 

526.8173 

Ar+l/2 

.72177 

526.5857 

_ 

.72213 

526.6007 

— 

Ar+1 

.72177 

526.2447 

— 

.72213 

526.2596 

- 

.72249 

526.2745 

526.2743 

The  optimum  value  of  a  changes  even  more  rapidly  in  the  fluorine  atom  going 
from  0.73732  for  F°  to  .72991  for  F_1.  Since  the  total  energy  depends  linear- 
Ly  on  a,  this  parameter  must  be  chosen  carefully. 

The  intersphere  exchange  coefficient  is  chosen  to  be  a  weighted  average 
of  the  atomic  exchange  parameters  from  the  two  constituents.  At  small  inter- 
nuclear  separations,  the  optimum  radius  for  an  atomic  sphere  frequently  places 
significant  amounts  of  charge  outside  that  atomic  sphere  -  charge  that  is 
still  strongly  associated  with  its  original  center  rather  than  being  trans¬ 
ferred  to  the  other  center  or  associated  with  the  molecular  binding  region. 

To  best  account  for  these  cases  the  weighting  coefficients  are  chosen  to 
reflect  the  origin  of  the  charge  in  the  intersphere  (or  outersphere 
region) , 


\  \  -  \  > +  \\  -  v> 

yintersphere  <  .  Q  o,  +  (  _  D) 

S1  1  S2  ^ 


where  (^s.  -  Q.°)  is  the  charge  lost  from  sphere  i  relative  to  its  atomic 

value  (or  ionic  value)  Q^°  and  as^  is  the  atomic  exchange  parameter  for  sphere 

i.  This  value  for  a.  ^  .  is  calculated  dynamically  -  it  is  updated 

intersphere  J  7  r 

after  each  iteration  in  the  self-consistent  calculation. 

While  for  heavy  atoms,  these  changes  in  the  exchange  parameter  would  be 
small,  the  a' s  for  small  atoms  vary  rapidly  with  2  (and  with  ionicity).  The 
correct  choice  of  the  exchange  parameters  influences  not  only  the  total  energy 
calculated  for  the  molecule  but  also  in  some  cases  affects  the  distribution  of 
charge  between  the  atomic  spheres  and  the  intersphere  region. 


B.  TRANSITION  PROBABILITIES 


The  electronic  and  vibrational-rotational  wavefunctions  of  a  pair  of 

states  can  be  used  to  calculate  transition  probabilities.  If  two  molecular 

states  are  separated  in  energy  by  an  amount  AE  =  hcv  (h  =  Planck's  constant, 

nm 

c  =  velocity  of  light,  V  =  frequency  in  wave  numbers),  the  semi-classical 
theory  of  radiation  (References  90  and  91)  yields  for  the  probability  of  a 
spontaneous  transition  from  an  upper  state  n  to  a  lower  state  m 

a  -  nm  Snm 

nm‘  3  t,4c3  -g7  (26) 


Here  A  is  the  Einstein  coefficient  for  spontaneous  transition  from  level  n  + 
nm 

m,  gn  is  the  total  degeneracy  factor  for  the  upper  state 


9n  =  (2-  8<r,A<)(2S,+  l)(2j'+l) 


(27) 


and  Snm  is  the  total  strength  of  a  component  line  in  a  specific  state  of 
polarization  and  propagated  in  a  fixed  direction.  A  related  quantity  is  the 
mean  radiative  lifetime  of  state  n  defined  by 


~  =  I 

Tn 


1  nm 


m<n 


(28) 


the  summation  being  over  all  lower  levels  which  offer  allowed  connections. 
The  intensity  of  the  emitted  radiation  is 


Inm-  A  E  nm  Nn  Anm 


(29) 


where  Nn  is  the  number  density  in  the  upper  state  n.  This  analysis 

assumes  that  all  degenerate  states  at  the  same  level  n  are  equally  populated, 

which  will  be  true  for  isotropic  excitation.  The  total  line  strength  S  can 

nm 

be  written  as  the  square  of  the  transition  moment  summed  over  all  degenerate 
components  of  the  molecular  states  n  and  m: 
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where  j  and  1  refer  to  all  quantum  numbers  associated  collectively  with 
upper  and  lower  electronic  states,  respectively. 

In  the  Born-Oppenheimer  approximation,  assuming  the  separability  of  elec 
tronic  and  nuclear  motion,  the  wavefunction  for  a  diatomic  molecule  can  be 
written  as 


where  (r,  R)  is  an  electronic  wavefunction  for  state  i  at  fixed  inter- 
nuclear  separation  R,  i(iv(R)  is  a  vibrational  wvefunction  for  level  v  and 
^JnA^0,  X*  ^  refers  to  the  rotational  state  specified  by  electronic  angular 
momentum  A,  total  angular  momentum  J  and  magnetic  quantum  number  M.  The 
representation  is  in  a  coordinate  system  related  to  a  space-fixed  system  by 
the  Eulerian  angles  (0,  x,  <t>) .  The  transition  moment  can  be  written, 
using  the  wavefunction  given  by  Equation  (31),  as 

Mi'  =/Vv'j'A'M'  {  ~e  +  !in}^v'’j"A’,M"  dTedrvdrr  (32) 


The  subscripts  e,  v  and  r  refer  to  the  electronic,  vibrational  and  rota- 
tional  wavefunctions  and  M  and  M  are  the  electronic  and  nuclear  electric 
dipole  moments,  respectively.  Integration  over  the  electronic  wavefunction, 
in  the  Born-Oppenheimer  approximation,  causes  the  contribution  of  the  nuclear 
moment  Mn  to  vanish  for  i  *  j.  The  electronic  dipole  moment  can  be  written 
(References  91  and  92)  in  the  form 


■?eIk=_t?exkK* 


where  the  primed  coordinates  refer  to  the  space  fixed  system,  the  coordin¬ 
ates  r^  refer  to  a  molecule-fixed  system  and  2>(0,  x»  is  a  group  rotation 
tensor  whose  elements  are  the  direction  cosines  related  to  the  Eulerian  rota¬ 
tion  angles  (0,  x>  $) •  Using  bracket  notation,  Equations  (32)  and  (33)  can  be 
combined  to  yield  for  the  transition  moment 


jv'j'A'M' 

M  j  j  =  M  = 


<^v'  -Z  er^  i  [6,x,  4>)  J"A"M'^> 


The  matrix  elements  <J'  A1  M*  |  (9,  x»  <t>)|j"  A"  M">  determine  the  group  se¬ 
lection  rules  for  an  allowed  transition  and  have  been  evaluated  for  many  types 
of  transitions  (References  36-38).  Summing  Equation  (34)  over  the  degenerate 
magnetic  quantum  numbers  M'  and  M"  we  have  from  Equation  (30) 


nw'J'A'  _  ^J'A'  n v' 
>>nm-  5mv"j"A"  '  «*j"A"  Pmv" 


J'  A’  . 


where  j»a»  is  the  Honl-London  factor  (References  39  and  40) 


is  the  band  strength  for  the  transition.  Combining  Equations  (27),  (29)  and 
(35),  we  have  for  the  intensity  of  a  single  emitting  line  from  upper  level 


Ta  C  nv  J  r_nv'j'A 

T  nv'j'  _4_..  lAEmv"j*J  smvV'A" 

1  n™  =  1  mv"j*  =  3  Nj'  •>  ~7~,  ‘ 

tl  C  W  r»(2  J  +  I  ) 


where  Nj  is  the  number  density  in  the  upper  rotational  state  J  and 

“n  =  <2-60iA)  (2S*  +  1)  is  the  electronic  degeneracy.  Taking  an 
’  n  v '  J 1 

average  value  of  vhj«  ^or  the  whole  band.  Equation  (37)  can  be  summed  to 
yield  the  total  intensity  in  the  (v1 ,  v")  band: 


/  nv'  l4  nv' 

'  .nv'j  _  _4_  /  Pmv" 

v"  -  3  Nv  •h4c3cor 


y 

where  N^,  *  N^,  is  the  total  number  density  in  the  upper  vibrational  level 

v'  and  where  we  make  use  of  the  group  summation  property 

I  =(2J,+  0  (39) 

J 

Comparing  Equations  (29)  and  (38),  we  have  for  the  Einstein  spontaneous 
transition  coefficient  of  the  band  (v1 ,  v") 


[A_  nv'  l3  nv' 
.  nv'  _  4  lAE  mv*J  Pmv 
Amv«=  3  t>*c  3u>n 


Similarly,  the  lifetime  of  an  upper  vibrational  level  v'  of  state  n  can  be 


written 


1  v  r*  A  nv 
Tn  =  L  4,  Amv  " 
n  m<n  v* 


where  the  summation  runs  over  all  v"  for  each  lower  state  m.  Equation  (40) 
can  be  cast  in  the  computational  form 

(2141759x10*)  rA1_nV'  /  ,1 3  nv'  (4; 

A mv "  ( sec  )-  —  [AEmv«(au)J  p|J,v//(au) 


where  AE  „  and  p  „  are  in  atomic  units.  It  is  also  often  convenient  to 
m  v  m  v 

relate  the  transition  probability  to  the  number  of  dispersion  electrons  needed 
to  explain  the  emission  strength  classically.  This  number,  the  f-number  or 
oscillator  strength  for  emission,  is  given  by 


fnm,v'v#  = 


me3  h2 


2e 


[ae"v'"]2 


A  nV' 

^mv" 


(43) 


The  inverse  process  of  absorption  is  related  to  the  above  development 
through  the  Einstein  B  coefficient.  Corresponding  to  Equation  (29),  we  have 
for  a  single  line  in  absorption 


I  nvA 

mV 


l°Ax 


J  K  [v)  d  v  -  AE  mn  B mn 

line  (vVj'V) 


(44) 


where  K(v)  is  the  absorption  coefficient  of  a  beam  of  photons  of  frequency 
v  and 


Q  dbv'j'A' 

Bmn=  Brny'^'A' 
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(45) 


is  the  Einstein  absorption  coefficient  for  a  single  line.  Summing  over  all 
lines  in  the  band  (v",  v' ) ,  assuming  an  average  band  frequency,  we  obtain 


where  Nv,t  =  j„  N^,,  is  the  total  number  density  in  the  lower  vibrational  state 
v" .  Corresponding  to  Equations  (42)  and  (43)  we  can  define  an  f-number  or 
oscillator  strength  for  absorption  as 
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In  computational  form,  Equation  (47)  becomes 
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where  AE  ^  and  v„  are  m  atomic  units.  Combining  Equations  (40)  and 

(43)  and  comparing  with  Equation  (47),  we  see  that  the  absorption  and  emission 

f-numbers  are  related  by 


f  mn,v"v' 


f  nm,v'  v" 


(49) 


Some  caution  must  be  observed  in  the  use  of  f-numbers  given  either  by  Equa¬ 
tion  (43)  or  (47)  since  both  band  f-numbers  and  system  f-numbers  are  defined 
in  the  literature.  The  confusion  arises  from  the  several  possible  band  aver¬ 
aging  schemes  that  can  be  identified. 

An  integrated  absorption  coefficient  (density  corrected)  can  be  defined 
from  Equation  (46)  as 
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(50) 


where  the  exponential  factor  corrects  for  stimulated  emission.  Equation  50 
can  be  written  in  terms  of  the  absorption  f-number  as 


Sv-y = 
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(51) 


Using  h  c/k  =  1.43880  cm-K° ,  we  obtain  a  computational  formula  for  the 

-2 .  -1 

integrated  absorption  coefficient  as  Svit  yi  (cm  ’atm  )  - 


2.3795 


1.43880 


exp 


V',  V 


(cm-1) 


T  (K‘ 


mn . 


(52) 


The  total  integrated  absorption  is  found  from 

c  -V  T  c (53) 

bT0TAL  ■  L  L  =»V  ,V 

v"  v' 

where,  under  normal  temperature  conditions,  only  the  first  few  fundamentals 
and  overtones  contribute  to  the  summations. 

The  developments  given  above  are  rigorous  for  band  systems  where  an  aver 
age  band  frequency  can  be  meaningfully  defined.  Further  approximations,  how¬ 
ever,  are  often  made.  For  example,  the  electronic  component  of  the  dipole 
transition  moment  can  be  defined  as 


This  quantity  is  often  a  slowly  varying  function  of  R  and  an  average  value 
can  sometimes  be  chosen.  Equation  (36)  can  then  be  written  approximately  in 
factored  form  as 


'f’mv"  “  qv'V.S  I  ^ji 

*»J 


where  q  ,  „  the  square  of  the  vibrational  overlap  integral,  is  called  the 
v  v  , 

Franck-Condon  factor,  eva^uatec^  at  some  mean  value  of  the  internuclear 

separation  R.  In  addition,  it  is  sometimes  possible  to  account  for  a  weak  R- 
dependence  in  M  by  a  Taylor  series  expansion  of  this  quantity  about  some 
reference  value,  usually  referred  to  the  (0,  0)  band.  We  have 


l  +  a(R-Ra/3)  +  b(R-Ra£)2+  ] 


Substituting  into  Equation  (56)  and  integrating  yields 


-Pmv*  -  qv'V.S  \  ®i'^  [  l  +  Q(Rv'v"-  Ra/3)  +  b(Rv'v" -  Ra£)2+  ]  |  (57) 


where 


(Rv'  v"  "  Ra/3)  = 


<V'|(R-Ra0)lv'> 


V|v*> 


is  the  R-centroid  for  the  transition  and 


'c 


is  the  R2-centroid.  Note  that  this  last  term  differs  (to  second  order)  from 
the  square  of  the  R-centroid.  An  alternate  procedure  can  be  developed  by 
evaluating  Equation  (54)  at  each  R-centroid,  Rvtvii.  Then 

•Pmv"  *  qvV'£>ii(Rv'v")|2  (60) 

i.J 


Equation  (60)  assumes  that  the  vibrational  wavefunction  product 
behaves  like  a  delta  function  upon  integration. 


v/V  ^v"  =  S(R—  Rv'v")  <y'|v"> 


The  range  of  validity  of  Equation  (60)  is  therefore  questionable,  particu¬ 
larly  for  band  systems  with  bad  overlap  conditions  such  as  oxygen  Schumann- 
Runge.  The  range  of  validity  of  the  R-centroid  approximation  has  been  exam¬ 
ined  by  Frazer  (Reference  41). 

The  final  step  in  calculating  transition  probabilities  is  the  determina¬ 
tion  of  Rj^(R),  the  electronic  dipole  transition  moment,  for  the  entire  range 
of  internuclear  separations,  R,  reached  in  the  vibrational  levels  to  be  con¬ 
sidered.  This  can  be  expressed  in  terras  of  the  expansion  of  Equation  (4) 


^j.(R)=Z  c!*c|,  <k  (R)|  Me|^(Rj) 

JJ  Is  > 


where  cyJ  and  c^  are  coefficients  for  and  respectively. 

An  analysis  similar  to  that  yielding  Equation  (10)  and  (11)  gives 


£ €P  <^MS I  ^ s P I <n(  lr k,R) | MeP|  n  Vrvk(rh,Rj> 


The  spatial  integral  in  Equation  (63)  reduces  to  one-electron  integrals 
equivalent  to  overlap  integrals,  and  the  evaluation  of  Equation  (63)  can  be 
carried  out  by  the  same  computer  programs  used  for  Equation  (11).  Programs 
for  evaluating  >^.^(R)  in  Equation  (62)  have  been  developed  at  UTRC  and  exam¬ 
ples  of  their  application  have  appeared  in  the  literature  (Reference  8). 

For  perturbed  electronic  systems,  the  transition  dipole  moment  will  have 
a  strong  R-dependence  and  R-centroid  or  other  apporximat ions  will  be  invalid. 
A  direct  evaluation  of  Equation  (36)  would  therefore  be  required  using  the 
fully-coupled  system  of  electronic  and  vibrational  wavefunctions  to  properly 
account  for  the  source  of  the  band  perturbations. 
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SECTION  3 


DISCUSSION  OF  RELATIVISTIC  METHODS 

For  heavy  atoms  (Z  >  30),  and  molecular  systems  built  from  heavy  atoms, 

relativistic  effects  become  increasingly  important  and  should  be  taken  into 

account  in  the  calculation  of  the  radial  wavefunctions.  The  implementation  of 

relativistic  effects  into  atomic  and  molecular  computer  codes  is  only  fairly 

recent  owing  to  the  increased  complexities  introduced  in  the  self-consistent 

field  (SCF)  procedure  and  the  greatly  increased  computer  time  required  for 

such  calculations.  Compared  with  the  non-relativistic  case,  the  Dirac- 

Hartree-Fock  (DHF)  method  requires  that  two  radial  functions,  G  ...  corre- 

nlj 

sponding  to  the  large  component  and  F  . .  corresponding  to  the  small  component 

nlj, 

must  be  calculated  for  each  of  the  two  possible  j  values.  Thus,  the  numerical 
work  of  a  DHF  relativistic  treatment  is  increased  by  nearly  a  factor  of  four 
over  the  nonrelativistic  case,  exclusive  of  increased  complexities  in  evalua¬ 
tion  of  the  terms  of  the  Hamiltonian.  In  view  of  this,  methods  that  have  been 
developed  to  date  for  molecular  systems  have  involved  the  use  of  model  poten¬ 
tials  to  represent  relativistic  effects. 

In  the  calculation  of  the  internal  energy  of  a  molecular  system  comprised 
of  n  electrons  and  N  nuclei,  and  considering  only  electrostatic  interactions 
between  the  particles,  we  have  for  the  total  Hamiltonian 
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where  m^,  m^,  M  ,  are  the  masses  of  the  electron,  atom  a  and  combined  system 
mass,  respectively.  Now  since  the  ratios  m  /m  and  m  />!_  are  both  small, 


(2x10  -  5x10  )  we  can  effect  a  separation  of  the  electronic  and  nuclear 

coordinates  treating  the  total  wavefunction  as  a  product  of  a  nuclear  and  an 
electronic  part.  We  have 


5n>  = 


where  ^(rn,  R^)  is  an  electronic  wavefunction  parametric  in  the  nuclear 
coordinates  as  given  in  Equation  (66)  and  are  nuclear  motion  wavefunc- 

tions  which  satisfy  (neglecting  terms  of  the  order  of  m&/ ma) 
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The  cross  term  in  •  7^  can  be  eliminated  by  a  proper  change  of  varia¬ 
bles  and  Equation  (67)  then  reduces  to  a  3N-3  dimensional  Schrodinger 
equation. 

For  most  systems,  where  the  velocity  of  motion  of  the  nuclei  is  slow 
relative  to  the  electron  velocity,  this  decoupling  of  electronic  and  nuclear 
motion  is  valid  and  is  referred  to  as  the  adiabatic  approximation.  Equation 
(66)  thus  defines  an  electronic  eigenstate  ,  R^),  parametric  in  the 

nuclear  coordinates,  and  a  corresponding  eigenvalue  E^(RN)  which  is  taken  to 
represent  the  potential  energy  curve  or  surface  corresponding  to  state  k. 

In  the  usual  jib  initio  method  for  calculating  the  electronic  properties 
of  a  molecular  system,  one  starts  from  a  zero-order  Hamiltonian  that  is  exact 
except  for  relativistic  and  magnetic  effects,  and  which  involves  the  evalua¬ 
tion  of  electronic  energies  and  other  relevant  quantities  for  wavefunct ions 
that  are  properly  antisymmetrized  in  the  coordinates  of  all  the  electrons. 

For  a  system  containing  n  electrons  and  M  nuclei,  the  zero-order  Hamiltonian 
depends  parametrically  on  the  nuclear  positions  and  is  of  the  form 


in2nM  z  i  M  z,  z,  | 
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where  z.  and  R.  are  the  charge  and  position  of  nucleus  i,  r^  is  the  position 
of  electron  j,  and  V.^  is  the  Laplacian  operator  for  electron  j.  All  quanti¬ 
ties  are  in  atomic  units,  i.e.  lengths  in  bohrs,  energies  in  hartrees  (1 
hartree  =  2  Rydbergs) . 

In  addition  to  the  electrostatic  contribution,  Je  ,  the  complete  Hamilton¬ 
ian  should  contain  additional  terms  which  correct  for  magnetic  interactions 
and  relativisitic  effects.  These  correction  terms  may  be  of  importance  in 
several  applications.  These  include: 

(1)  calculation  of  the  probability  of  making  a  transition  from  one  quan¬ 
tum  state  to  another  in  high-momentum  collisions  such  as  those  that 
can  occur  in  hot  atom  or  heavy  atom  chemical  dynamics  experi¬ 
ments; 

(2)  determination  of  the  interaction  energy  in  heavy  nuclei  systems  such 
as  Cs2  and  U0+,  which  exhibit  open-shell  structure  on  both  nuclei  at 
infinite  internuclear  separations; 

(3)  calculation  of  the  intermolecular  forces  between  free  radicals, 
electronically  excited  states  of  molecules  with  open-shell  struc¬ 
ture,  and  long  molecular  conformations  of  possible  biological 
interest . 

A.  B RE IT- PAULI  HAMILTONIAN 

The  relativistic  correction  terms  to  the  usual  electrostatic  Hamiltonian 

have  been  derived  through  order  az ,  where  a  is  the  fine  structure  constant, 

and  are  often  referred  to  as  the  Breit-Pauli  (Reference  42)  Hamiltonian  terms. 

This  Hamiltonian  has  been  derived  by  Bethe  and  Salpeter  (Reference  43)  for  a 

two-electron  system  and  has  been  generalized  to  the  many-electron  system  by 

Hirschfelder ,  et  al  (Reference  44)  and  Itoh  (Reference  45).  In  the  absence  of 

external  electric  or  magnetic  fields  we  can  represent  these  correction  terms 

as  follows.  Let  s.  and  p.  =4-4.  denote  the  operators  for  the  spin  and  linear 
J  J  1  J 
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moment  of  electron  j,  respectively.  Then  the  generalized  Breit-Pauli  Hamil¬ 
tonian,  correct  to  terms  of  Q( et2/M),  can  be  written  as: 


*^BP"  +a^LL  +  ^ SS  +  ^LS  +  + 


where  ^e  is  given  by  Equation  (68)  and  the  correction  terms  can  be 
expressed  as  follows: 
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The  first  correction  term L,  represents  the  magnetic  orbit-orbit  coupling 
terms  of  the  electrons  arising  from  the  interaction  of  the  magnetic  fields 
created  by  their  motion.  The  second  term,^  gives  the  spin-spin  magnetic 

O  b 

coupling  terms  which  are  often  quite  appreciable.  For  r.  =  0,  only  the 

jk 

delta-function  contribution  survives  which  represents  the  Fermi-contact  spin 
interaction.  The  third  ,  is  usually  the  largest  in  magnitude  and 

Lb 

represents  the  spin-orbit  interaction  between  the  spin  and  magnetic  moment  of 
each  electron  and  the  spin-other  orbit  interaction,  which  represents  the 
coupling  of  the  spin  of  one  electron  with  the  magnetic  moment  of  a  different 
electron.  The  termJtyj(p,  corrects  for  variation  of  the  electron  mass  with 
velocity  and  the  represents  electron  spin  terms  identified  by 

Dirac  which  appear  to  have  no  classical  analogue. 


Aside  from  the  spin-orbit  terra,^^,  usually  only  the  last  term,^//^, 
(often  called  the  Darwin  correction  term)  andjf^,  the  mass-velocity  term, 
are  retained  in  the  Hamiltonian,  yielding  the  so-called  Pauli  approximation 
(Reference  43) . 

The  eigenfunctions  of  the  Hamiltonian  represented  by  Equation  (69)  are 
four-component  Dirac  spinors  which  may  be  expressed  as: 


/nkm  ~  y 


Pnk(r)  Xkm(0,<*>) 

iQnk(r)  X-km<G,*) 


where  4>)  are  products  of  spherical  harmonics  and  Pauli  spinors  and 

Pn^(r),  Qnk(r)  represent,  respectively,  the  large  and  small  components  of  the 
radial  wave  equation.  The  exact  solution  of  the  Breit-Pauli  Hamiltonian  has 
only  been  given  for  one-  and  two-electron  atomic  systems  (Reference  43)  owing 
to  the  complexity  of  the  operators  for  the  general  n-electron  case.  For  a 
molecular  system,  Kolos  and  Wolniewicz  (Reference  46)  have  calculated  the 
relativistic  corrections  to  H2  using  Equation  (69)  to^ta2).  No  heavier  mole¬ 
cular  systems  have  been  treated  using  the  full  Breit-Pauli  Hamiltonian. 
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B.  APPROXIMATE  TREATMENTS 


Although  the  Breit-Pauli  HamiLtonian  given  in  Equation  (69)  can  formally 
be  employed  in  a  molecular  system,  both  the  multiplicity  of  terms  and  the 
difficulty  of  evaluation  of  the  resultant  molecular  integrals  has  precluded 
its  general  use  to  date.  For  atomic  systems,  various  approximate  methods  of 
solution,  within  a  Hartree-Fock  or  multiconfiguration  Hartree-Fock  framework, 
have  been  proposed  for  atoms  (References  47-51).  In  most  of  these  methods,  a 
restricted  Hamiltonian  which  includes  only  the  one-electron  Dirac  terms  is 
usually  employed.  The  contributions  of  the  Breit  operators  for  spin-magnetic 
interactions  and  velocity  retardation  are  then  calculated  as  first-order  per¬ 
turbations  using  the  zeroth-order  Dirac  relativistic  wavefunct ions . 

An  even  more  approximate  method  for  incorporating  the  major  relativistic 
effects  has  been  proposed  by  Cowan  and  Griffin  (Reference  52).  In  this 
method,  the  mass-velocity  and  Darwin  terms,  written  in  terms  of 

the  Pauli  equation  fof  one-electron  atoms,  are  simply  added  to  the  usual  non- 


relativistic  Hamiltonian  operator.  In  addition,  the  spin-orbit  terms, , 

Li  O 


are  omitted,  thereby  reducing  the  system  of  equations  to  a  single  form  repre¬ 


senting  the  description  of  the  major  component  wavefunction,  P^^r),  evaluated 


at  the  center-of-gravity  of  the  spin-orbit  states.  The  rationale  for  this 
approximation  lies  in  the  observation  that  detailed  atomic  calculations  using 
the  complete  DHF  method  have  indicated  that,  even  for  an  atom  as  heavy  as 
uranium,  less  than  1  percent  of  the  total  charge  is  described  by  the  small 
component  radial  wavefunct ions . 

The  resulting  equations  have  the  form: 


_di  +  MVJJ  +  vi(r)  +  (r)  +  HDj  (r) 

dr2  r 2 


G  nl  =  e'nf  G'n£(r) 


(76) 


where  the  mass-velocity  and  Darwin  terms  take  the  form 
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and  a  -  1/137.036  is  the  fine  structure  constant.  The  spin-orbit  term  is 
omitted  in  Equation  (76)  and  thus  these  equations  represent  center  of  gravity 
radial  functions  averaged  over  the  two  possible  total  angular  momentum  quantum 
numbers.  Equation  (76)  represents  (apart  from  the  neglect  of  spin-orbit 
effects)  the  relativistic  corrections  to  first  order  in  a2.  A  more  accurate 


analysis  of  heavy  atom  energy  levels  and  spectra  is  available  through  the  use 
of  the  radial  functions,  G^^tr),  found  from  Equation  (76),  and  a  first-order 
perturbation  calculation.  Cowan  and  Griffin  (Reference  52)  have  illustrated 
the  utility  and  accuracy  of  such  an  approach. 

Recently,  Wood  and  Boring  (Reference  53)  have  adapted  this  approximate 
relativistic  method  to  the  local  exchange  problem  and  have  implemented  the 
so'ution  of  Equation  (76)  within  the  context  of  the  multiple  scattering  Xa 
raetnod  (Reference  22).  The  central  field  Hamiltonian  is  modified  to  include 
mass-velocity  and  Darwin  terras,  given  by  Equations  (77)  and  (78),  in  the 
sphere  surrounding  each  atomic  center.  The  intersphere  region  in  the  multiple 
scattering  approach  (constant  potential  region)  is  treated  nonrelativistically 
since  charge  in  this  region  is  far  from  a  nucleus  and  is  screened  by  the 
charge  concentrated  around  the  atomic  centers.  The  matching  conditions  for 
continuity  of  the  wavefunction  at  the  sphere  boundaries  permits  any  necessary 
charge  transfer  between  the  relativistic  intra-atomic  regions  and  the  nonrela- 
tivistic  interatomic  constant  potential  regions.  For  an  atom,  the  Wood-Boring 
treatment  reduces  to  the  Dirac-Slater  local  exchange  method,  but  with  the 
neglect  of  spin-orbit  terms. 


The  implementation  of  Equation  (76)  into  existing  nonrelat ivist ic  multi¬ 
ple  scattering  molecular  codes  is  facilitated  by  a  change  in  the  dependent 
variable,  G^^Cr),  to  eliminate  the  first  derivative  of  the  wavefunct ion , 
illustrated  in  Equation  (78).  The  usual  Numerov  method  of  solution  can  then 
be  applied  to  the  central  field  problem;  the  only  new  requirement  being  the 
numerical  tabulation  of  the  first  and  second  derivatives  of  the  potential  at 
each  grid  point  in  the  integrations.  These  derivatives  are  computed  only  once 
for  each  complete  SCF  cycle  and  thus  the  total  required  computer  time  for  a 
typical  problem  is  not  significantly  increased  as  compared  with  a  nonrelativ- 
istic  calculation.  A  complete  self-consistent  program  incorporating  this 
method  has  been  developed  at  UTRC.  Our  code  has  been  tested  by  repeating 
calculations  for  the  U  and  Pu  atoms  (Reference  52),  where  we  find  excellent 
agreement  with  the  more  exact,  but  cumbersome,  Dirac-Slater  calculations . 
Results  for  molecular  calculations  have  recently  been  reported  by  Boring  and 
Wood  for  UF6  and  U02++  (References  54  and  55).  These  calculations  were  car¬ 
ried  out  to  illustrate  the  shifts  in  the  valence  levels  for  such  systems  re¬ 
sulting  from  relativistic  effects.  The  total  energy  was  not  of  principal 
concern. 

We  have  recently  reported  (Reference  56)  the  first  all-electron  calcula¬ 
tion  of  the  potential  energy  curves  for  a  molecule  (Hg2+)  built  from  atoms 
which  exhibit  significant  relativistic  effects.  This  study  illustrated  that 
reliable  total  energies  are  obtainable  through  a  relativistic  multiple  scat¬ 
tering  density  functional  treatment,  provided  care  is  taken  to  optimize  poten¬ 
tial  match  and  overlap  criteria  for  such  systems.  This  study  formed  the  basis 
of  the  computational  scheme  that  we  have  employed  here  for  the  uranium/oxygen 
system. 

C.  EFFECTIVE  CORE  MODELS 

It  is  well  known  from  chemical  experience  that  the  outermost  valence 
electrons  contribute  most  to  determining  the  chemical  properties,  especially 
spectroscopic  properties,  of  molecules.  The  core  electrons  remain  essentially 
unchanged  from  their  atomic  form  except  for  internuclear  separations  of  the 
order  of  the  charge  radii  of  the  outer  core  region  or  less,  wherein  core  po¬ 
larization  effects  may  become  important.  Since  the  computational  time 


required  for  ab  initio  calculations  of  electronic  structure  goes  up  at  least 
quadrat ical ly  with  the  number  of  electrons  in  the  system,  there  have  been  many 
attempts  to  replace  the  more  tightly  bound  core  electrons  with  simple  one- 
electron  effective  potentials  (References  57-65).  Concurrent  with  elimination 
of  an  explicit  treatment  of  the  core  electrons,  a  transformation  of  the 
valence  orbital  basis  is  required  to  insure  that  the  lowest  valence  orbital  of 
each  symmetry  has  a  nodeless  radial  form,  since  it  is  well  known  that  the 
lowest  energy  eigenfunction  for  a  local  potential  must  be  nodeless  (Reference 
66). 

Typical  of  the  several  effective  core  models  that  have  been  reported  is 
that  due  to  Kahn,  et  al  (Reference  64)  whereby  an  effective  core  potential  is 
described  in  terms  of  angularly  dependent  projection  operators  as 

U c°re  =  U  (r)  +  I  I  |lm>[ufcore  (r)  -ULcore(r)]  <im|  (79) 

I  m 

where  L  is  taken  at  least  as  large  as  the  highest  angular  momentum  orbital 
occupied  in  the  core.  The  term  U  core(r)  represents  the  effective  Coulomb  and 
exchange  potential  felt  by  the  valence  electrons.  The  second  term  essentially 
accounts  for  the  repulsive  potential  between  valence  and  core  electrons  for 
each  symmetry  i.  The  only  non-local  character  exhibited  by  a  potential  of  the 
form  of  Equation  (79)  arises  from  the  ^.-dependence  which  can  be  cast  in  terms 
of  one-electron  integrals  between  the  core  and  valence  orbitals.  Explicit 
two-electron  terms  connecting  core  and  valence  orbitals  are  thus  avoided  which 
greatly  simplifies  the  calculation  of  matrix  elements  of  the  effective  Hamil¬ 
tonian.  The  potential  given  by  Equation  (79)  can  be  compared  with  the  gener¬ 
alized  Phi l 1 ips-Kleinman  pseudo-potential  (Reference  59). 


(80) 


UCOrC  =  I  (2JC  -  kc)  +  VC0 


where  Jc  and  Kc  represent  the  core  orbital  Coulomb  and  exchange  operators 


.CO 


and  V  is  a  complicated  non-local  operator  which  guarantees  core-valence 
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orthogonal Lty .  Since  the  P-K  core  orbitals  must  simultaneously  be  eigenfunc- 

CO 

tions  of  both  the  core  and  valence  Hartree-Fock  Hamiltonians,  V  ,  in  general, 
contains  complicated  two-electron  terms  and  limits  the  usefulness  of  Equation 
(80)  over  a  full  ab  initio  treatment. 

The  prescription  of  Kahn  can  be  implemented  by  analytically  fitting  a 
nodeless  pseudo-orbital,  xn£»  to  a  linear  combination  of  numerical  or  analytic 
Hartree-Fock  orbitals  determined  from  a  full  self-consistent  treatment  of  the 
core  electrons,  using,  for  example,  the  multiconfiguration  Hartree-Fock  code 


of  Froese-Fischer  (Reference  67).  The  components,  U 


(r),  of  Equation 


(79)  are  then  defined  implicitly  from  the  Schrodinger  equation 


whereby 


2 

[  2  F“  +  Ui  (r)  +  2JvQ)-KvQl]  Xni_€nl  Xn<£ 


2  J vO  I  + 


Kval]  XnJ 


Equation  (81)  can  be  extended  to  relativistic  systems  in  several  ways. 

Kahn,  et  al  (Reference  68)  suggest  an  approximate  treatment  of  adding  only 

the  mass-velocity  and  Darwin  terras  to  the  usual  electrostatic  Hamiltonian  and 

to  determine  approximate  HF  orbitals  in  the  manner  prescribed  by  Cowan  and 

Griffin  (Reference  52).  Equation  (81)  is  then  used  to  determine  an  effective 
c  or© 

U  (r)  such  that  e  are  the  eigenvalues  of  the  CG  approximate  relativis- 

x.  n* 

tic  solution  and  ^  are  curve-fitted  to  the  CG-HF  orbitals.  In  this  treat¬ 
ment,  the  xn£  represent  approximate  solutions  to  the  major  component  wavefunc- 
tion,  Pn£,  determined  at  the  center-of-gravity  of  the  spin-orbit  states. 

Lee,  et  all  (Reference  65)  adopt  a  somewhat  more  complicated  treatment  in 
which  the  spin-orbit  operator  is  added  to  the  usual  electrostatic  Hamiltonian, 


in  addition  to  the  mass-velocity  and  Darwin  terms  retained  in  the  Cowan- 
Griffin  treatment.  The  large  component  eigenfunctions  of  a  full  Dirac- 
Hartree-Fock  treatment  of  the  atom,  as  given,  for  example,  by  Desclaux  (Refer¬ 
ence  69)  are  then  curve-fitted  in  a  manner  similar  to  the  Kahn  treatment  but 
include  the  additional  index  for  the  particular  spin-orbit  state,  Xn£j •  Use 
of  these  eigenfunctions  in  a  molecular  system  fits  more  naturally  into  a  (J-J) 
coupling  scheme  whereas  the  xn^  determined  using  the  Kahn  method  are  more 
easily  represented  using  A-s  coupling. 

Although  these  effective  core  models  can  often  accurately  describe  an 
atomic  eigenvalue  sequence,  including  even  high-lying  electronically  excited 
states  (Reference  70),  there  are  inherent  difficulties  in  their  application  to 
molecular  environments,  where  the  maximum  angular  momentum  component  of  the 
valence  shell  orbitals  may  often  exceed  the  highest  Jl-value  component  retained 
in  Equation  (82).  This  is  particularly  true  for  valence  orbitals  which  exhi¬ 
bit  strong  changes  from  atomic  form  through  hybridization  with  higher  angular 
momentum  orbitals  or  through  the  addition  of  more  compact  polarizaiton  terms. 
In  either  case,  since  the  relativistic  terms  are  now  all  buried  in  a  fixed 
rather  than  a  dynamic  relativistic  operator,  only  static  core  effects  are 
imposed  in  determining  the  shape  of  the  valence  molecular  orbitals.  Relativ¬ 
istic  effects  between  valence  electrons  and  shielding  effects  of  the  core  by 
the  valence  electrons  are  therefore  neglected  in  these  effective  core  treat¬ 
ments.  In  addition,  the  models  obviously  break  down  completely  when  the 
nuclei  are  brought  together  to  dimensions  such  that  core  overlap  and  polariza¬ 
tion  effects  become  significant.  Unfortunately,  calculations  to  date  seem  to 
indicate  that  such  effects  begin  to  set  in  for  internuclear  separations  of  the 
order  of  equilibrium  bond  lengths. 


SECTION  4 


DISCUSSION  OF  RESULTS 

A  preliminary  study  of  the  electronic  structure  of  U,  U+,  UO  and  U0+  had 
previously  been  undertaken  by  UTRC.  Both  Hartree-Fock-Slater  (Reference  71) 
and  screened  hydrogenic  (Reference  72)  wavefunct ions  had  been  used  to  find 
optimum  7s,  5f,  and  6d  orbitals.  Only  small  differences  in  the  orbitals  were 
found  between  these  two  approaches  and  the  calculated  orbitals  were  fitted  to 
LCAO  expansions  using  an  effective  nuclear  charge  Hamiltonian.  The  energies 
of  the  various  possible  sublevels  of  the  ground  state  configuration  for  U 
(7s  5f  6d;  158  terms)  and  U  (7s  5f  ;  17  terms)  were  then  calculated  using 

these  effective  nuclear  charges.  These  calculated  term  levels  are  in  good 
agreement  with  the  experimental  studies  of  Kiess,  et  al .  (Reference  73)  for  U 
and  the  group-theoretical  treatment  of  Elliott,  et  al .  (Reference  74)  for  U+. 
Relativistic  corrections  and  spin-orbit  splitting  factors  were  estimated  using 
perturbation  theory  based  on  a  one-electron  model.  It  was  found  that  L-S 
coupling  is  a  reasonable  approximation  for  many  multiplets  of  these  systems, 
provided  there  is  little  configuration  mixing  or  promotion  from  the  7s 
orbitals,  which  are  found  to  have  strong  relativistic  corrections.  Departure 
from  pure  L-S  coupling  among  the  5f  electrons  was  estimated  to  be  less  than  10 
percent  (Reference  75).  A  complete  (J-J) coupled  analysis  of  U  and  U+  was  not 
carried  out  in  these  preliminary  studies. 

An  analysis  of  the  electronic  structure  of  UO  and  U0+  using  a  relativis¬ 
tic  formulation  was  undertaken  by  UTRC  under  this  Contract.  The  possible  low- 
lying  molecular  states  are  shown  in  Tables  1  and  5-  Preliminary  calculations 
were  performed  for  several  states  of  UO  and  U0+  and  for  the  ground  state  of 
UO2 ,  U0j+  and  U02++.  A  brief  summary  of  the  results  of  these  molecular  calcu¬ 
lations  follows. 

UO 

Electronic  structure  calculations  were  carried  out  for  this  system  using 
a  relativistic  density  functional  formalism.  Only  a  selected  group  of 
symmetries  was  studied.  Our  calculations  indicate  that  the  lowest  symmetry  of 
UO  is  derived  from  the  (A,S)  coupled  ^1  state  and  has  the  following  principal 
molecular  orbital  occupancy: 
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'  '  3 

The  J  3ir  16  1 4>  group  derives  from  the  5f  atomic  configuration  in  the  U 


atom  and  is  quartet  coupled.  We  have  found  that  the  5I  state  of  UO  is  the 
ground  state  but  that  several  other  symmetries,  including  5K  and  3I,  are  low- 
lying.  The  results  are  shown  in  Figure  1  which  indicates  that  only  states  of 
triplet  or  quintet  multiplicity  are  bound  for  this  system.  An  examination  of 
the  structure  of  these  low-lying  states  of  UO  indicates  a  near  total  charge 
transfer  to  U+20“2  for  short  (equilibrium)  internuclear  separations.  Thus 
only  the  molecular  states  in  the  lower  group  shown  in  Table  1  are  likely  to  be 
bound.  This  would  yield  42  bound  molecular  states  arising  from  u[5l]  +  o[3p] 
and  39  repulsive  states.  Since  (J-J)  coupling  is  surely  a  better 
approximation  for  UO,  these  two  manifolds  of  states  will  be  optically 
connected  and  many  pre-dissociation  paths  of  the  type: 


UO  +  hv  *  UO  >  U  +  0 


(84) 


. 

are  possible.  Here  UO  is  a  vibrat ional ly  excited  low-multiplet  state  of  UO 

**  ...  .  . 
and  UO  is  a  dissociating  state.  The  predicted  optical  absorption  should  be 


strong  since  the  transfer  is  from  highly  ionic  states  to  neutral  valence 
states  of  UO.  Since  UO  has  a  large  dissociation  energy  (7.87  eV),  both  one 
photon  and  two  photon  solar  excitation  processes  are  possible. 

A  vibrational  analysis  of  the  ft  =  5  ground  state  of  UO  was  carried  out 
using  a  Hulbert-Hirshfe lder  (Reference  76)  fit  to  our  calculated  potential 
curves.  This  fit  yields  an  equilibrium  internuclear  distance  of  1.90  A  and  a 
fundamental  vibrational  constant  (we)  of  859  cm-1.  The  spin-orbit  interac¬ 
tion  was  calculated  using  U+2  atomic  splitting  parameters.  No  explicit  two- 
center  effects  are  included.  Our  calculated  spectroscopic  data  are  compared 
in  Table  2  with  the  work  of  Krauss  and  Stevens  (Reference  77)  and  estimates 
based  on  experimental  data  for  similar  systems.  The  agreement  is  well  within 
the  uncertainty  of  the  calculations  or  experimental  estimates. 
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An  analysis  of  the  LWIR  emission  from  UO  was  carried  out  based  on  the 
ground  electronic  state.  These  calculations  should  also  be  representative 
of  the  LWIR  emission  from  other  low-lying  electronic  states  since  they  exhibit 
similar  ionicities.  Our  calculated  f-numbers,  including  fundamentals  and 
overtones,  are  given  in  Table  3.  These  data  (v*  >  v")  are  given  for  the 
lowest  30  vibrational  levels.  We  have  also  included  the  emission  wavelengths 
for  each  transition.  For  UO,  our  calculated  f-number  for  the  1-0  transition 
is  4.86  x  10~5  at  X  =  12.04  y.  In  Table  4,  we  present  the  calculated  inte¬ 
grated  band  absorption  coefficients  as  a  function  of  temperature.  All  over¬ 
tone  contributions  have  been  included  in  these  band  absorption  coeffi¬ 
cients  . 

UO+ 

Detailed  searches  of  several  symmetries  of  UO+  were  carried  out  to  deter¬ 
mine  the  ground  molecular  state  of  this  system.  Our  calculations  indicate 
that  the  lowest  symmetry  of  U0+  is  derived  from  the  (A,  S)  coupled  4I  state 
and  has  the  following  principal  molecular  orbital  occupancy: 

4I  [la2  2 a2  3 a2  4 a2  lx4  2tt4  !  3*  1<J>  16  I]  (85) 

I  I 

Again  we  find  that  the  3x  l<f>  16  group  is  quartet  coupled  in  the  ground 
state  but  a  second  manifold  of  states  for  the  U0+  system,  which  exhibit 
doublet  coupling  of  these  electrons,  lies  about  2  eV  above  the  ground  state. 

An  apparent  gap  in  the  density  of  states  for  UO+  is  found  between  these  two 
groups . 

Calculations  for  U0+  proved  to  be  much  more  complex  than  UO  owing  to  the 
presence  of  at  least  two  low-lying  dissociation  limits  of  U+  +  0.  As  indi¬ 
cated  in  Table  5,  the  ionic  U++0“  structures  will  mix  in  all  multiplicities 
with  molecular  ion  states  arising  from  U+[4l]  +  0[^p],  U+[^l]  +  o[^p],  and 
U  [  k]  +  o[  p].  Thus,  this  system  may  not  exhibit  the  repulsive  valence 
states  that  were  found  for  UO. 

A  vibrational  analysis  of  the  Q  =  9/2  ground  state  of  UO"1-  was  carried  out 
using  a  Hulbert-Hirschf elder  (Reference  76)  fit  to  our  calculated  potential 


I  III 


t  '  . ’  splittings  were  derived  from  atomic  parameters  for  the 
•  1  Is  ,m  e<|ii  i  1  i  hr  i  um  intern  uclear  distance  of  1.84  A  and  a 
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t  md  in.  f t  i  i  it  i  n  il  mist  ant  of  890  cm'*.  These  data  are  compared  in 

Tab’..'  witt  ftwi  h;  dated  estimates,  since  there  are  no  experimental  data 
avai  lab’..-.  The  ig  t  .■•■m.-nt  ber  ween  our  work  and  that  of  Krauss  and  Stevens 
(Reference  ’ 7 ‘  is  bos  satisfactory  than  in  the  case  of  UO  but  still  well 
within  the  uncertainty  ot  the  several  c a  1 c u  1  at  i  ons . 

\  perturbative  treatment  for  calculating  the  density  of  states  in  uranium 
molecules  is  available  through  the  use  of  ligand  field  theory.  The  basic 
concept  relies  on  the  assumption  that  the  structure  and  density  of  the  5f 
electrons  in  uranium  and  its  ions  remain  unchanged  in  a  molecular  environment. 
Our  calculated  multiplet  splitting  for  the  5f  electrons  in  U+  is  given  in 
Table  6.  ’’"his  clearly  indicates  that  the  4I  component  of  5^  is  low-lying  and 
that  the  lowest  doublet  manifold  lies  ~  2.2  eV  above  the  ground  state.  If 
this  splitting  carries  over  to  the  U0+  ion  without  much  change,  the  following 
solar  pumped  process  is  possible: 


UO* [4I j  +  hv  +  uo**[2h] 
e  +  U0+*(2H]  ♦  U  [5L]  +  0[3P] 


(86) 


Dissociative-recombination  is  not  energetically  possible  from  the  ground  4I 
state  of  U0+  but  it  is  possible  from  the  excited  2H  and  higher  states.  The 
optical  connections  between  these  two  manifolds  of  UO*  states  will  be  studied 
in  future  work. 

An  analysis  of  the  emission  characteristics  for  the  ground  state  of  U0+ 
indicates  an  oscillator  strength  for  emission  C f  ^ 0 )  of  5.17  x  10-5  at  A  = 

11.3  p.  A  complete  analysis  of  our  calculated  LWIR  emission  for  U0+  is  given 
in  Table  7  for  the  lowest  30  vibrational  levels.  In  Table  8,  we  present  our 
calculated  integrated  band  absorption  coefficients  for  U0+.  A  detailed 
analysis  of  the  vibrational  -  rotational  spectra  and  band  absorption/emission 
characteristics  for  these  systems,  UO  and  U0+ ,  is  currently  in  progress. 

Our  calculated  LWIR  emission  for  U0+  is  typical  of  that  for  a  highly 
ionic  r-'Jtal  oxide.  We  predict  strong  emission  from  the  fundamentals  of  U0+  in 


the  wavelength  region  11.3  ~  14  y.  Since  this  system  exhibits  weak 
anharmonic i ty ,  we  fine  the  overtones  down  in  intensity  by  several  orders  of 
magnitude.  However,  the  first  excited  state  of  U0+  (4H)  lies  at  ~  2800  cm-1 
in  our  calculations  with  a  predicted  electronic  oscillator  strength  of  ~  1  x 
10-5  for  the  4I  -  4H  transition.  An  analysis  of  the  spin-orbit  splitting  of 
the  excited  4H  manifold  has  not  yet  been  carried  out  but  the  predicted  result 
is  that  the  lowest  4H7/2  component  should  lie  at  ~  1500  cm-1  with  a  similar 
oscillator  strength  (~  1  x  10-5)  for  the  5Ig/2  ~  ^7/2  transition.  Thus  the 
electronic  and  vibrational  manifolds  for  UO+  are  highly  overlapped  above  the 
second  vibrational  level  of  the  ground  5Ig/2  state.  Since  the  density  of 
electronic  states  of  U0+  is  very  large  above  ~  2.2  eV,  we  predict  that  strong 
solar  pumping,  followed  by  intense  radiation  in  the  region  0.6  <  X  <  11.3  y 
should  occur  for  this  system.  This  conclusion  is  similar  to  that  reached  by 
Krauss  and  Stevens  (Reference  77)  based  on  their  MCSCF  analysis  of  the  UO+ 
system.  Since  the  excited  electronic  states  of  U0+  lying  in  the  region  of 
strong  solar  flux  ( .4  -  .7  y)  exhibit  shifted  equilibrium  internuclear 
separation  from  that  of  the  ground  4I  state,  we  predict  efficient  conversion 
of  solar  photons  to  IR  photons  for  this  system. 

uo2 

Our  calculations  of  the  triatomic  species,  U02 ,  have  indicated  two  low- 
lying  electronic  states  in  symmetry.  Their  principle  molecular  orbital 
occupancies  are: 


JHU  [  1  ag^  lou2  2og2  litu^  2ou2  3dg2  2ttu4  lug4  3ou^  i  16g  l<i>u] 
5IU  [log2  lau2  2og2  lwu^  2ou2  3ag2  2ttu^  lug^  3ou  l  3uu  16g  1  cpu ] 


4  ,..4  ->„..2  ! 


3  =  +4  =  5 

The  lowest  H  state  corresponds  to  an  ionicity  of  0  U  0  ,  whereas  the  I 
u  u 

state  indicates  back-bonding  from  the  oxygen  ligands  into  the  stable  5^  sub¬ 
structure  of  uranium.  A  third  low-lying  5I  structure  is  also  found  which 

8 

corresponds  to  a  30^  4o^  promotion  out  of  the  second  configuration  given 

above.  More  detailed  calculations  of  this  system  are  indicated  but  our 


preliminary  data  indicate  that  U02  should  exhibit  an  integrated  band 
absorption/emission  coefficient  in  the  LWIR  of  between  1.5  -  2.0  times  that 
calculated  for  UO. 

uo2+ 

A  preliminary  analysis  of  UO  +  has  been  carried  out  in  D  symmetry. 

I  2 

McGlynn  and  Smith  (Reference  78)  have  proposed  a  <5^  ground  state  for  this  ion 
based  on  simple  molecular  orbital  arguments.  Their  set  of  one-electron 
orbital  energies  is  based  on  a  maximum  overlap  criterion  that  is  empirical  in 
character.  More  modern  calculations  of  the  actinide  series  atoms  suggest  that 
a  U+5[0_2]2  or  U+3[0_1]  2  structure  should  be  the  most  stable  configuration. 

In  terms  of  MO's,  the  lowest  predicted  electronic  states  (D^)  would  be: 


log2  lou2  2ag2  liru^  2cru2  3ag2  2ttu^  lirg^  3au2 


1  <i>u 


or 


(88) 


log2  lau2  2ag2  lnu^  2ou2  3og2  2itu^  lirg^  3au  16g  1  <pu  : 


Further  calculations  should  be  carried  out  to  determine  which  of  these 

configurations  has  lower  energy  and  to  examine  several  other  possible  low- 

.  .  .  2 
lying  excited  states.  Our  calculations  to  date  indicate  that  the  symmetry 

is  the  ground  state.  If  the  highly  ionic  2<J>u  state,  as  predicted  by  McGlynn, 
is  in  fact  the  ground  state,  the  radiation/absorption  characteristics  of  U02 
should  be  approximately  twice  as  strong  as  U0+.  No  real  quantitative  informa¬ 
tion  about  the  absorption  spectrum  of  UO  +  is  available  at  present  but  McGlynn 

.  2  2  2  ~ 
predicts  a  strong  ^-j/2 u  +  ^5/2g  transiti°n  at  X  -  0.4y.  Again,  this  transi¬ 

tion  lies  in  the  peak  solar  flux  region. 

There  appear  to  be  many  strong  optical  absorption  bands  in  this  system 
corresponding  to  molecular  orbital  transitions  between  MO's  derived  from  f+d 
uranium  atom  symmetries.  There  are  also  some  large  ligand  to  central  uranium 
atom  charge  transfer  effects  corresponding  to  transitions  between  the  two 
types  of  configurations  listed  above.  Our  preliminary  data  indicate  that 
U02+  should  possess  a  rich  LWIR  spectra  and  should  exhibit  many  strong  elec¬ 
tronic  bands  (f  ~  10”2  -  10-1)  beginning  at  ~  500  nm. 
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An  extensive  set  of  calculations  for  the  ground  state  of  UO^  was 
carried  out  in  symmetry.  The  ground  electronic  state  has  the  dominant 

molecular  orbital  configuration: 

[log2  lou2  2og2  lnu^  2ou2  3og2  2ttu^  lirg^  3ou2]  (89) 

We  find  an  equilibrium  internuclear  separation  of  2.06  A,  a  value  somewhat 
larger  than  that  corresponding  to  the  ion  in  solution  (~  1.7  A).  At  present 
it  is  uncertain  whether  this  represents  a  deficiency  in  the  calculations  or 
whether  significant  back-bonding  of  electron  charge  to  the  central  uranium 
atom  is  occurring  in  the  gas  phase  UC>2++  species.  This  would  result  in  a 
lengthing  of  the  U-0  bond  and  a  corresponding  decrease  in  the  U-0  bond 
strength. 

The  first  excited  electronic  state  of  UO  ++  corresponds  to  a  3ou  +  1  g 
electron  promotion.  This  lies  at  ~  3.6  eV  and  we  therefore  predict  very  weak 
solar  pumping  for  this  species.  The  LWIR  analysis  should  be  carried  out, 
however,  since  the  highly  ionic  nature  of  UC>2++  may  give  rise  to  a  strong 
absorption/emission  character  for  the  asymmetric  stretch  vibrational  mode. 


SECTION  5 


RECOMMENDATIONS 

The  low-lying  molecular  states  of  UO+  arising  from  ^1  (5f^7s2), 
^L(5f^7s6d),  and  ^K(5f^7s6d)  of  U+  are  given  in  Table  2.  Additional  low-lying 
states  of  U0+  arise  from  the  ^M(5f^6d7p)  state  of  U+  and  from  the  doubly 

n  n  /  "X 

ionized  5L(5f  6d)  and  5 1 ( 5 f  7s)  and  triply  ionized  I ( 5 f  )  states  of  the 
uranium  ion.  Further  calculations  of  the  term  manifold  for  this  molecular  ion 
are  needed.  In  our  preliminary  studies,  a  separate  (higher  lying)  chemistry 
was  found  for  U+  when  the  5f  electrons  were  doublet  spin-coupled. 

Verification  of  this  structure,  using  a  relativistic  framework  for  this  ion, 
is  an  area  of  proposed  research. 

Finally,  theoretical  studies  of  U0++ ,  UC^  and  the  uranyl  (UC>2++)  ion 
should  be  undertaken.  The  UO^  ion  has  the  indicated  ground  state  configura- 
t  ion: 

2$u  [log2  lau2  2ag2  liru^  2ou2  3ag2  2itu^  lug^  3au2  Hu]  (90) 

where  the  composition  of  the  various  molecular  orbitals  is  given  below: 


Molecular  Orbital 


Composition 


log 

lou ,  2og 
Iitu,  2ou 

3og,  2ttu,  3au,  litg 
l<t>u 


6s(U) 

2s(0) 

6p(U) 

Mainly  2p(0)  +  5f(U)  +  6d(U) 
5f(U) 


The  ordering  and  energetics  of  the  low-lying  symmetries  of  U02+  is 
uncertain.  A  detailed  study  of  this  ion,  examining  both  relativistic  and 
spin-coupling  effects  is  indicated  as  an  area  of  proposed  research.  Studies 
in  both  linear  and  bent  conformations  should  be  carried  out.  These  studies 
should  include  an  analysis  of  the  IR  radiation  properties  of  the  ground 
electronic  state  and  any  possible  low-lying  electronic  transitions  in  the  IR 
region . 

For  most  of  these  proposed  calculations  a  density  functional  relativistic 
framework  could  be  employed.  Spin-orbit  splittings  and  the  degree  of  mixing 
from  (A-S)  to  (J-J)  coupling  could  be  calculated  using  perturbation  tech¬ 
niques.  Our  final  relativistic  wavefunct ions  could  then  be  used  to  calculate 


both  the  vibration-rotation  and  electronic  transitions  in  these  ions,  that 
fall  in  the  LWIR  region  of  interest.  The  low-lying  electronic  states  up  to 
about  3  eV  excitation  should  be  included  in  this  study. 

The  focus  of  these  studies  should  continue  to  be  on  the  U+,  U0+  and  U02  + 
systems  but  a  limited  study  of  U0++  and  U02++  is  indicated.  An  examination  of 
the  thermodynamics  of  the  uranium/oxygen  system  as  given  in  Figure  2,  shows 
that  U+  and  U0+  are  terminal  ions  and  therefore  potentially  important 
radiators  at  high  altitudes,  whereas  U02+  is  more  important  at  altitudes  of 
significant  02  concentration.  The  importance  of  U0++  and  U02++  is  not  well 
known  owing  to  the  uncertainty  in  the  thermodynamics  of  these  ions.  It  is 
anticipated  that  continued  theoretical  studies  will  clear  up  many  of  the 
current  uncertainties  in  the  radiation  charac terist ics  of  these  systems. 

Follow-on  studies  should  also  include  an  analysis  of  the  several  possible 
kinetic  pathways  that  have  been  suggested  by  our  work  to  date.  These  include 
photoexcitation  and  photodissociation  of  UO  as: 


U+0~  +  hv  ♦  UO* 
UO*  +  hv  U  +  0 


(91) 


where  UO*  represents  the  excited  neutral  valence  states  of  UO.  Other  poten¬ 
tially  important  kinetic  studies  include  dissociative-recombination  routes 
from 

e  +  [U0+]*  +  U*  +  0  (92) 


where  [UO  ]  is  an  electronically  excited  state  of  U0+  which  exhibits  doublet 
coupling  for  the  5f  electrons  and  charge  transfer  processes  such  as 


In  addition,  the  dissociative-recombination  routes 
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Figure  2.  Thermodynamic  stability  of  known  uranium/oxygen  species. 


Table  1.  Molecular  states  of  UO. 


Separated  Atom 


Molecular  States 


U  +  0 


5L  (SfVfid)  +  3p.(2p4) 


E«>  *  0.000  eV 
(81  states) 


3  3 

£"(1),  V(2),  V(  1),  V(2),  V(l), 

7E+(2),  3n(3),  5n(3)t  7H(3) ,  3a(3),  5ac 
7A(3),  3*(3),  5*(3),  7*(3) ,  3r(3),  5r(3) 
7r(3),  3H(3) ,  5H(3) ,  7H(3),  3I (3) ,  51(3) 
71(3),  3K (3) ,  5K(3),  7K(3),  3A(2) ,  5A(2) 
7A(2),  3M(1),  5M(1),  7M(i) 


+  3Pg(2p*)  VU),  V«) .  *,-(2).  »!♦„),  V(2).  7: 


Eos  -  0.096  eV 
(72  states) 


3n(3),  5n(3),  7n(3),  3A(3),  5A(3),  7A(3), 
3*(3),  5*(3),  7»(3),  3F(3) ,  5F(3),  7r(3) , 
3H(3),  5H(3),  7H(3),  3I (3) ,  51(3),  ?I(3), 
3K(2),  5k(2),  7K(2),  3A(1)  ,  5A(l),  7A(1) 


Table  1. 

Molecular  states  of  UO.  (Continued) 

Separated  Atom 

Molecular  States 

U+  +  0" 

4Iu(5f37s2)  +  2Pu(2p5) 

3Z+(1),  V(l),  V(2),  V(2),  3 

(42  states) 

5n(3) ,  3A (3) ,  5A(3),  30 ( 3 ) »  5$(3) 

3r(3),  5r(3),  3H(3) ,  5H(3),  3I (2) 

5I(2) ,  3K(1),  5k(i) 

U4"4  +  0= 

5Lu(5f36d)  +  1Sg(2P6) 

5Z~(1),  5n(i) ,  5a ( i) ,  54>(1),  5r(i 

(9  states) 

5H(l),  5I (1) ,  5K(1),  5A(1) 

Table  2.  Calculated  spectroscopic  data  for  UO/UO+ 


Up J  5i  ] 

U0+  [  4 1  1 

coe  (cm_l) : 

820 

(matrix) 

949 

(SCF) 

836 

(scaled  from  ThO) 

935 

(MCSCF) 

863 

(SCF)* 

925 

(MCSCF:  Q  =  9/ 

845 

(MCSCF) 

890 

(RDF-UTRC) 

859 

(RDF-UTRC) 

u)eXe(cm~l) : 

2.7 

(H-H  potential) 

2.7 

(H-H  potential) 

12.0 

(3l  excited  state) 

Re(A): 

1.84 

(matrix) 

1.83 

(SCF) 

1.88 

(SCF) 

1.84 

(MCSCF) 

1.84 

(MCSCF-Krauss) 

1.84 

(RDF-UTRC) 

1.90 

(RDF-UTRC) 

*SCF  and  MCSCF 

results 

,  see  Reference  77. 
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Table  3-1.  Calculated  oscillator  strengths  (fviv»)  for  the  vibrational-rotational  transition  of  UO(X5I-x5l) 


oscillator  strength 
wavelength,  microns 


Table  3-2.  Calculated  oscillator  strengths  (fv«v")  for  the  vibrational-rotational  transition  of  UO(X^I  X^I) 


oscillator  strength 
wavelength,  micron! 


oscillator  strength 


ulatcd  oscillator  strengths  (lv*v")  for  the  vibrational-rotational  trans 


oscillator  strength 
wavelength,  microns 


oscillator  strength 
wavelength,  microns 


Table  3-6.  Calculated  oscillator  strengths  (fv'v")  f°r  the  vibrational-rotational  transition  of  UO(x5l-x5l) 


oscillator  strength 
wavelength,  microns 


Table  4. 

Total  integrated 

absorption  coefficients  for 

UO  ground  state 

Absorption  Coefficient, 

-2  -1 
cm  atm 

Temperature,  °K 

Fundamental 

First 

Overtone 

Second 

Overtone 

Third 

Overtone 

(Wavelength,  y) 

(12.04) 

(6.04) 

100. 

1156.665 

3.728 

273.15 

1156.581 

3.823 

300. 

1156.538 

3.869 

500. 

1155.947 

4.479 

1000. 

1153.399 

6.974 

1500. 

1150.375 

9.893 

2000. 

1147.191 

12.951 

2500. 

1143.911 

16.082 

3000. 

1140.551 

19.265 

4000. 

1133.599 

25.747 

5000. 

1125.879 

32.279 

Table  5 

Separated  Atom 
lf*~  +  0 

4I  (5f37a2)  +  3PR(2p4) 

U  O 

£.  -  6.11  eV 
(63  states) 


^Of^d)  +  3P  (2p4) 

Em  »  6.146  eV 
(81  states) 


*Ku(5f37s6d)  ♦  3Pg(2p4) 


Molecular  states  of  U0+. 


Molecular  States 


V(l),  V(2).  V(l).  V(2),  V(l). 

6z*( 2).  2n(3),  4n(3),  6n(3),  2a(3),  4ai 


6M3), 

,  26(3) , 

4*(3), 

6*(3), 

,  2r(3), 

4r  (3: 

6r(3), 

,  ^O), 

4H(3), 

6H(3), 

2K2), 

4i(2: 

6I(2), 

,  ^(1), 

4K(1), 

6K(1) 

4  _ 
1(1) 

,  V<2) 

1 ,  V(l)  ,  V 

(2),  V(l), 

4n'3), 

6H(3). 

sn(3) , 

4A(3) , 

6A(3), 

®A(3) 

4f<3), 

6*(3) , 

**(3), 

4r(3), 

6r(3), 

®r(3) 

4H(3), 

®H(3), 

®H(3), 

41(3) , 

6I(3) , 

®1(3) 

4K(3), 

VS). 

®K(3>, 

4A(2), 

6A(2), 

®A(2) 

^(l). 

^1(1). 

®M(1) 

V<2) 

.  4r+(i),  6r(2),  V 

(l).  ®r 

"(2), 

4H(3), 

6n(3), 

®n(3) , 

4A(3) , 

®A(3), 

®A(3) 

4*(3). 

*♦(3). 

®*(3) , 

4r(3), 

6r(3), 

®r(3) 

4H(3), 

*H(3), 

®H(3), 

4I(3). 

6I(3), 

®I(3) 

4K(2), 

®K<2), 

®K(2) , 

4A(1), 

6A(1), 

®A(1) 

Em  -  6.223  eV 
(72  states) 


Table  5.  Molecular  states  of  U0+  (Continued) 


Separated  Atom 
U++  +  Q- 

5Lu(5f36d)  +  2Pu(2p5) 
(54  states) 


U4~f+  +  0~ 

4Iu(5f3)  +  1Sg(2p6) 
(8  states) 


Molecular  States 


4Z+(1),  4E“(2) ,  6Z+(1),  6z“(2) ,  4n 
6n(3),  4A (3) ,  6A(3),  4$(3),  6$(3), 
4r (3) ,  6r(3),  4H(3),  6h(3),  4I(3), 
6I(3),  4K(3),  6K(3),  4A (2) ,  6A(2), 
^(1),  ^(1) 


4z"(i) ,  4n(i) ,  4a(i) ,  4*(i),  4r(i) 

4H(l),  41(1),  4K(1) 


Calculated  oscillator  strengths  (fv'v")  for  the  vibrational-rotational  transition  of  UO+CA^I-A^I) 
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Table  7-5.  Calculated  oscillator  strengths 


h 


i-H 

00 

vO 

co 

o 

rH 

CO 

o 

m 

m 

• 

CM 

>3* 

m 

CO 

NO 

CM 

NO 

CM 

rH 

00 

• 

CM 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

ON 

rH 

CO 

NO 

rH 

<r 

CO 

co 

CM 

rH 

CM 

co 

rH 

m 

NO 

00 

ON 

rH 

CO 

CM 

o 

o 

o 

o 

rH 

l 

m 

rH 

rH 

nO 

ON 

o 

CM 

NO 

00 

NO 

m 

NO 

<r 

CO 

On 

CO 

CO 

CM 

CO 

CM 

ON 

m 

00 

rH 

r- 

00 

o 

NO 

CM 

NO 

rH 

NO 

CO 

CM 

CM 

On 

CM 

CM 

rH 

rH 

rH 

vD 

00 

ON 

rH 

CO 

CO 

CO 

o 

o 

o 

rH 

rH 

rH 

rH 

rH 

ON 

ON 

CO 

NO 

in 

CO 

CM 

CO 

H 

CM 

i— 1 

rH 

co 

m 

CO 

rH 

rH 

00 

On 

NO 

m 

rH 

m 

co 

CM 

CM 

r- 

CM 

CM 

rH 

m 

rH 

rH 

00 

ON 

rH 

CO 

CO 

co 

<r 

o 

o 

rH 

rH 

rH 

rH 

rH 

rH 

00 

uo 

ON 

co 

o 

CO 

<r 

rH 

ON 

CO 

CO 

o 

CM 

rH 

r* 

in 

rH 

CM 

00 

NO 

00 

ON 

CO 

<r 

CO 

rH 

CM 

m 

CM 

CM 

rH 

rH 

CM 

rH 

ON 

rH 

on 

rH 

CO 

CO 

<1* 

o 

rH 

rH 

rH 

rH 

rH 

rH 

m 

CO 

NO 

rH 

00 

CM 

NO 

O 

ON 

CM 

O 

ON 

rH 

00 

CO 

IO 

rH 

rH 

NO 

00 

rH 

NO 

rH 

rH 

CM 

<r 

rH 

H 

CM 

CM 

rH 

CO 

rH 

rH 

rH 

rH 

CO 

rH 

o 

r*. 

rH 

CM 

ON 

00  rH 

CO 

00 

CO 

r>. 

m 

00 

rH 

00 

00 

CM 

m 

o  <r 

vC 

CM 

m 

rH 

o 

CM 

CM 

ON 

rH 

rH 

.H 

ON  rH 

rH 

m 

rH 

rH 

rH 

00 

ON 

rH 

r*- 

in 

CM 

m 

• 

• 

• 

• 

CO 

rH 

rH 

m 

ON 

NO 

CO 

NO 

00 

ON 

CM 

CM 

m 

rH 

o 

ON 

CO 

rH 

CM 

rH 

rH 

rH 

co 

o 

rH 

NO 

CO 

ON 

NO 

00 

m 

CM 

m 

CO 

rH  CM 

rH 

rH 

oo 

rH 

NO  »H 

rH 

00 

NO 

NO 

00 

oo 

rH 

m 

o 

00 

ON 

CM 

ON 

CM 

rH 

CM 

o 

rH 

O 

rH 

ON 

m 

in 

CM 

rH 

m 

CM 

CO 

CO 

CM 

o 

rH 

oo 

o 

CM 

rH 

rH 

rH 

rH 

rH 

CO 

rH 

o 

On 

rH 

oo 

m 

ON 

O 

ON 

CM 

00 

rH 

CO 

in 

CM 

00 

rH 

Tt'fV 


18-16  2.83-16  4.25-16  4.27-15  1.54-16  2.63-15  2.31-15  1.68-14  1.36 

67  0.71  0.76  0.81  0.87  0.93  1.01  1.10  1.20 


Table  7-6.  Calculated  oscillator  strengths  (fv«v")  for  the  vibrational-rotational  transition  of  U0+ 


*  oscillator  strength 
**  wavelength,  microns 


Table  8.  Total  integrated  absorption  coefficients  for  U0+  ground  state 


Absorption  Coefficient,  cm  atm 


Temperature,  °K 

Fundamental 

First 

Overtone 

Second 

Overtone 

Third 

Overtone 

(Wavelength,  p) 

(11.30) 

(5.67) 

(3.79) 

(2.85) 

100. 

1230.827 

3.810 

.023 

.000 

273.15 

1230.764 

3.882 

.024 

.000 

300. 

1230.729 

3.921 

.024 

.000 

500. 

1230.215 

4.458 

.974 

.000 

1000. 

1227.852 

6.782 

.061 

.001 

1500. 

1224.990 

9.549 

.117 

.002 

2000. 

1221.961 

12.462 

.197 

.004 

2500. 

1218.825 

15.449 

.302 

.008 

3000. 

1215.540 

18.485 

.433 

.013 

4000. 

1208.159 

24.668 

.777 

.032 

Total 

1234.660 

1234.670 

1234.674 

1235.647 

1234.696 

1234.658 

1234.624 

1234.584 

1234.471 

1233.636 


5000 


1199.135 


30.928 


1231.355 
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ATTN:  Tech  Library 

Battel le  Memorial  Institute 
ATTN:  R.  Thatcher 
ATTN:  STOIAC 
ATTN:  H.  Lamuth 

BDM  Corp 

ATTN:  T.  Neighbors 
ATTN:  L.  Jacobs 

Berkeley  Rsch  Associates,  Inc 
ATTN:  C.  Prettie 
ATTN:  J.  Workman 
ATTN:  S.  Brecht 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Boeing  Aerospace  Co 

ATTN:  MS/87-63,  D.  Clauson 

Boeing  Co 

ATTN:  G.  Hall 
ATTN:  S.  Tashird 

Boston  College 

ATTN:  W.  Grieder 
ATTN:  E.  Hegblom 

Boston  College 

2  cy  ATTN:  Chairman,  Dept  of  Chemistry 
2  cy  ATTN:  Chairman,  Dept  of  Physics 

BR  Communications 

ATTN:  J.  McLaughlin 

University  of  California  at  Riverside 
ATTN:  J.  Pitts,  Jr 

University  of  California  at  San  Diego 
ATTN:  H.  Booker 

California  Insittute  of  Technology 
ATTN:  J.  Ajello 

Cal  span  Corp 

ATTN:  C.  Treanor 
ATTN:  J.  Grace 
ATTN:  M.  Dunn 
ATTN:  W.  Wurster 

Charles  Stark  Draper  Lab,  Inc 
ATTN:  A.  Tetewski 
ATTN:  J.  Gilmore 
ATTN:  D.  Cox 

University  of  Colorado 

ATTN:  C.  Lineberger-JILA 
ATTN:  G.  Lawrence-LASP 

Columbia  University 

ATTN:  Security  Officer  for  H.  Foley 

Computer  Sciences  Corp 

ATTN:  F.  Eisenbarth 

Comsat  Labs 

ATTN:  D.  Fang 
ATTN:  G.  Hyde 

Concord  Sciences 

ATTN:  E.  Sutton 

Cornell  University 

ATTN:  D.  Farley,  Jr 
ATTN:  M.  Kelly 

University  of  Denver 

ATTN:  Sec  Officer  for  D.  Murcray 

University  of  Denver 
ATTN:  E.  Van  Zyl 

E-Systems,  Inc 

ATTN:  R.  Berezdivin 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Electrospace  Systems,  Inc 
ATTN:  H.  Logston 
ATTN:  P.  Phillips 

Environmental  Rsch  Inst  of  Michigan 
ATTN:  1RIA  Library 

EOS  Technologies,  Inc 
ATTN:  B.  Gabbard 
ATTN:  W.  Lelevier 

General  Electric  Co 


ATTN: 

R.  Edsall 

ATTN: 

P.  Zavitsanos 

ATTN: 

C.  Zierdt 

ATTN: 

R.  Juner 

ATTN: 

A.  Steinmayer 

General  Electric  Co 

ATTN:  G.  Mi 11  man 

General  Research  Corp 
ATTN:  B.  Bennett 

GEO  Centers,  Inc 

ATTN:  E.  Marram 

GTE  Communications  Products  Corp 
ATTN:  R.  Steinhoff 

GTE  Communications  Products  Corp 
ATTN:  J.  Concordia 
ATTN:  I.  Kohl  berg 

Harris  Corp 

ATTN:  E.  Knick 

Honeywell,  Inc 

ATTN:  A.  Kearns,  MS924-3 
ATTN.  G.  Terry,  Avionics  Dept 

Horizons  Technology,  Inc 


ATTN 

R.  Kruger 

HSS,  Inc 

ATTN 

M.  Shuler 

ATTN 

D.  Hansen 

IBM  Corp 

ATTN 

H.  Ulander 

Institute  for  Defense  Analyses 

ATTN 

H.  Gates 

ATTN 

E.  Bauer 

ATTN 

J.  Aein 

ATTN 

H.  Wolfhard 

International  Tel  &  Telegraph  Corp 
ATTN:  Tech  Library 

International  Tel  &  Telegraph  Corp 
ATTN:  G.  Wetmore 

IRT  Corp 

ATTN:  H.  Mitchell 

JAYCOR 

ATTN:  J.  Sperling 
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DEPARTMMENT  OF  DEFENSE  CONTRACTORS  (Continued 


Johns  Hopkins  University 


ATTN 

T.  Evans 

ATTN 

K.  Potocki 

ATTN 

C.  Meng 

ATTN 

J.  Newland 

ATTN 

P.  Komiske 

ATTN 

J.  Phillips 

Kaman 

Sciences  Corp 

ATTN 

J.  Jordano 

Kaman 

Sciences  Corp 

ATTN 

E.  Conrad 

Kaman  Tmepo 

ATTN 

W.  Schuleter 

ATTN 

DAS I AC 

ATTN 

W.  McNamara 

ATTN 

B.  Gambill 

Kaman  Tempo 

ATTN 

DAS  I  AC 

Litton  Systems,  Inc 

ATTN:  B.  Zimmer 


Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  J.  Reagan 
ATTN:  J.  Cladis 
ATTN:  J.  Kumer 
ATTN:  M.  Walt 
ATTN:  R.  Sears 
ATTN:  J.  Perez 
ATTN:  B.  HcCormac 

Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  Dept  60-12 
2  cy  ATTN:  D.  Churchill 

MIT  Lincoln  Lab 

ATTN:  0.  Towle 
ATTN:  V.  Vitto 
ATTN:  N.  Doherty 

MA/COM  Linkabit  Inc 

ATTN :  I .  Jacobs 
ATTN:  A.  Viterbi 
ATTN:  H.  Van  Trees 

Magnavox  Govt  &  Indus  Electronics  Co 
ATTN:  G.  White 

Maxim  Technologies,  Inc 
ATTN:  J.  Marshall 
ATTN:  R.  Morganstern 
ATTN:  E.  Tsui 

McDonnell  Douglas  Corp 
ATTN:  R.  Halprin 
ATTN:  Technical  Library  Svcs 
ATTN:  W.  Olson 

Meteor  Communications  Corp 
ATTN:  R.  Leader 

Pacific-Sierra  Rsch  Corp 
ATTN:  E.  Field,  Jr 
ATTN:  F.  Thomas 
ATTN:  H.  Brode,  Chairman  SAGE 


Mission  Research  Corp 


ATTN 

R.  Dana 

ATTN 

R.  Bogusch 

ATTN 

D.  Archer 

ATTN 

F.  Guigliano 

ATTN 

M.  Scheibe 

ATTN 

R.  Hendrick 

ATTN 

G.  McCartor 

ATTN 

R.  Bigoni 

ATTN 

R.  Kilb 

ATTN 

D.  Knepp 

ATTN 

F.  Fajen 

ATTN 

S.  Gutsche 

ATTN 

D.  Sowle 

ATTN 

Tech  Library 

ATTN 

P.  Fischer 

4  cy  ATTN 

C.  Lauer 

Mitre  Corp 

ATTN 

A.  Kyrmel 

ATTN 

MS  J104,  M.  Dresp 

ATTN 

G.  Harding 

ATTN 

C.  Callahan 

Mitre  Corp 

ATTN 

W.  Hall 

ATTN 

W.  Foster 

ATTN 

J.  Wheeler 

ATTN 

M.  Horrocks 

Nichols  Research  Corp,  Inc 
ATTN:  R.  Burns 

Pennsyvania  State  University 

ATTN:  Ionospheric  Rsch  Lab 

Photometries,  Inc 

ATTN:  I.  Kofsky 

Physical  Dynamics,  Inc 


ATTN 

A. 

Thompson 

Physical  Dynamics,  Inc 

ATTN 

J. 

Secan 

ATTN 

E. 

Fremouw 

Physical  Research,  Inc 

ATTN 

J. 

Devore 

ATTN 

R. 

Deliberis 

ATTN 

T. 

Stephens 

Physical  Science  Lab 

ATTN 

W. 

Berning 

Physical  Sciences,  Inc 

ATTN 

R. 

Taylor 

ATTN 

D. 

Greem 

ATTN 

T. 

Rawlings 

ATTN 

G. 

Caledonia 

ATTN 

K. 

Wray 

Physics  International  Co 
ATTN:  Tech  Library 

University  of  the  Commonwealth,  Pittsburgh 
ATTN:  M.  Biondi 
ATTN:  F.  Kaufman 
ATTN:  W.  Fite 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Princeton  University 
ATTN:  Librarian 

R&D  Associates 

ATTN:  C.  Grei finger 
ATTN:  M.  Gantsweg 
ATTN:  R.  Turco 
ATTN:  F.  Gilmore 
ATTN:  H.  Ory 
ATTN:  R.  lindgren 
ATTN:  W.  Karzas 
ATTN:  W.  Wright 
ATTN:  G.  Stoyr 
ATTN:  P.  Haas 

R&D  Associates 

ATTN:  J.  Rosengren 
ATTN:  B.  Yoon 


Rand  Corp 
ATTN 
ATTN 
ATTN 


P.  Davis 
C.  Crain 
E.  Bedrozian 


Rand  Corp 

ATTN:  B.  Bennett 

Riverside  Rsch  Institute 
ATTN:  V.  Trapani 

Rockwell  International  Corp 
ATTN:  R.  Buckner 

Rockwell  International  Corp 
ATTN:  S.  Guilici 

Science  Applications,  Inc 
ATTN:  D.  Sachs 
ATTN:  C.  Smith 
ATTN:  D.  Hamlin 
ATTN:  L.  linson 
ATTN:  E.  Straker 

Science  Applications,  Inc 
ATTN:  E.  Hyman 
ATTN:  J.  Cockayne 

Science  Applications,  Inc 
ATTN:  M.  Cross 

Space  Data  Corp 

ATTN:  S.  Fisher 

Spectral  Sciences,  Inc 
ATTN:  F.  Bien 

SRI  International 

ATTN:  C.  Hulburt 

Toyon  Rsch  Corp 

ATTN:  J.  Ise 
ATTN:  J.  Garbarino 
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SRI  International 

ATTN:  G.  Smith 
ATTN:  J.  Casper 
ATTN:  M.  Baron 
ATTN:  V.  Gonzales 
ATTN:  W.  Jaye 
ATTN:  J.  Petri  ekes 
ATTP  A.  Burns 
ATTN:  C.  Rino 
ATTN:  R.  leadabrand 
ATTN:  J.  Vickrey 
ATTN:  G.  Price 
ATTN:  R.  Tsunoda 
ATTN:  R.  Livingston 
ATTN:  W.  Chesnut 
ATTN:  D.  McDaniels 
ATTN:  D.  Neil  son 

Stewart  Radiance  Lab 
ATTN:  J.  Ulwick 

Swerling,  Manasse  &  Smith,  Inc 
ATTN:  R.  Manasse 

Technology  International  Corp 
ATTN:  W.  Boquist 

Teledyne  Brown  Engineering 
ATTN:  F.  Leopard 
ATTN:  N.  Passino 
ATTN:  MS-12,  Technical  Library 

TRW  Electronics  &  Defense  Sector 
ATTN:  G.  Kirchner 
ATTN:  R.  Plebuch 

Utah  State  University 

Attention  Sec  Control  Ofc  for 
ATTN:  C.  Wyatt 
ATTN:  A.  Steed 

ATTN:  K.  Baker,  Dir  Atmos  &  Space  Sci 
ATTN:  D.  Burt 
ATTN:  D.  Baker 


L.  Jensen,  Elec  Eng  Dept 


Visidyne,  Inc 
ATTN: 
ATTN:  I 
ATTN: 
ATTN: 
ATTN:  I 
ATTN:  l 


J.  Carpenter 
W.  Reidy 
T.  Degges 
0.  Shepard 
H.  Smith 
C.  Humphrey 


Wayne  State  University 
ATTN:  R.  Kummler 

Wayne  State  University 
ATTN:  W.  Kauppila 

United  Technologies  Rsch  Ctr 
2  cy  ATTN:  H.  Michels 
2  cy  ATTN:  R.  Hobbs 
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